o 

nmM 

^  j 

ly  AFCRL-64-677 

c> 


/ 

c/ 

MC  63-Ba-Rl 

GAB  DYNAMICS  DF 
HIGH-ALTITUDE  ROCKET  PLUMES 

HENRY  L.  ALDCN 
ROQER  H.  HABERT 

LJNTRACT  NO.  Art9<62B)-32BD 
PROJECT  NO.  49SA  TASK  NO.  49B4D1 

FINAL  REPORT  v/ 

JULY  1964 


Q 

< 


PREPARED  FDR 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  STATES  AIR  FORCE  < 

BEDFORD.  MASSACHUSETTS 


h-iA- 

/a:7u. 


MITHRAS,  Inc. 

MICROWAVE  AND  THERMAL  RADIATION  BENSORS 


!" 


.  _  .  300  PUTNAM  AVCNUe.  CAMBRIOOf.  MASS.  02139 


1 


Requests  for  additional  copies  by  Agencies  of  the  Department 
of  Defense)  their  contractors)  or  other  /government  agencies 
should  be  directed  to  the 

Defense  Documentation  Center  (DDC) 

Cameron  Station 
Alexandria,  Virginia  22314 

Department  of  Defense  contractors  must  be  established  for  DDC 
services  or  have  their  "need-to-know"  certified  by  the  cognizant 
military  agency  of  their  project  or  contract. 

All  other  persons  and  drgahi^atidhs  should  apply  •'<  the 

Ui  Si  Department  of  Commerce 
Office  of  Techhical  Services 
Washington,  D.  C.  20230 


MITHRAS,  INC. 

380  Putnam  Avenue 
Cambridge,  Massachusetts 
02139 


AFCRL  64-677 


MC  63-80-Rl 

GAS  DYNAMICS  OF 
HIGH-ALTITUDE  ROCKET  PLUMES 

Henry  L.  Aldeh 
Roger  H.  Habert 


Contract  No.  AF19(628)-3280 
Project  No.  4984 
Task  No.  498401 
ARPA  No.  42 


FINAL  REPORT 
July  1964 


This  research  is  part  of  Project  DEFENDER  sponsored  by  the 
Advanced  Research  Projects  Agency,  Department  of  Defense. 


Prepared 

for 

Air  Force  Cambridge  Research  Laboratories 
Office  of  Aerospace  Research 
United  States  Air  Force 
Bedford,  Massachusetts 


FOREWORD 

This  report  was  prepared  by  MITHRASi  Inc.,  of  Cambridge, 
Massachusetts,  for  the  Upper  Atmosphere  Physics  Laboratory,  Air 
Force  Cambridge  Research  Laboratories,  Hanscom  Field,  Bedford, 
Massachusetts, under  Contract  AF19(6Z8)-3280.  The  work  was  initiated 
by  Dr.  N.  W.  Rosenberg  and  monitored  by  him  and  D'r.  D.  Golomb. 
The  research  was  supported  by  Advanced  Research  Projects  Agency 
as  part  of  Project  DEFENDER. 

The  investigations  whose  results  are  reported  were  conducted 
during  the  period  1  July  1963  to  30  June  1964.  They  were  directed  by 
Mr.  Jacques  Ai  F.  Hill  and  carried  but  by  Mr.  Henry  L.  Alden,  and 
Mr.  Roger  Hi  Habert.  This  final  report  was  written  by  Mn  Alden. 

The  work  reported  here  represents  a  continuation  of  a  program 
begun  during  the  previous  year  under  Contract  AFl 9(628)-415.  This 
report  concludes  the  work  bn  Contract  AF19(628)-3 280. 
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ABSTRACT 


This  report  describes  methods  for  calculation  of  the 
velocities  and  thermodynamic  properties  of  flow  fields 
through  and  around  high  altitude  rocket  plumes.  The  main 
effort  is  directed  at  obtaining  a  gotxi  representation  of  both 
the  jet  flow  from  rocket  motors  and  the  shock-mixing 
layers  in  the  frontal  regions;  the  calculation  then  proceeds 
downstream  using  standard  methodsi  The  work  may  be 
divided  in  three  parts:  the  development  of  methods  to  cal¬ 
culate  the  flow  fields  throughout  plumes  using  the  principles 
of  inviscid  continuum  mechanics;  the  study  of  the  shock¬ 
mixing  layer  structure  at  the  nose  of  the  plume  using  the 
general  Navier -Stokes  equations  (i.e. »  accounting  for  vis¬ 
cosity*  heat  conduction  and  diffusion);  and  an  estimate  Of 
validity  of  continuum  mechanics  for  describing  high  altitude 
plumes*  as  inferred  from  a  study  of  shock  wave  formation 
at  high  altitudes.  A  comparison  with  available  observations 
shows  good  agreement  with  overall  plume  dimensions  as 
determined  by  the  flow-field  analysis.  It  may  be  concluded 
that  the  analysis  provides  a  sound  gas -dynamic  description 
which  maybe  used  as  a  basis  for  investigation  of  the  plasma 
and  chemi- electric  properties  of  high  altitude  plumes. 
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1.  INTRODUCTION 


As  the  exhaust  gases  released  by  a  missile  in  powered  flight 
expand  into  the  atmosphere,  a  complex  of  fluid  dynamic,  chemical  and 
radiative  phenomena  occur.  These  are  associated  with  a  large  visible 
plume  which  forms  around  the  rocket  and  trails  behind  it,  as  shown 
by  Rosenberg  (1961).  The  radiative  phenomena  are  "observables"  and 
have  been  subjected  to  the  most  intensive  observation.  In  order  to 
understand  how  they  arise  and  what  form  they  may  take  in  new  situa¬ 
tions,  it  is  necessary  to  study  thi  chernistry,  plasma  properties  and 
fluid  dynamics  Of  the  plumes. 

The  work  reported  here  is  Concerned  with  the  fluid  dynamic 
aspects  of  missile  tra'ls,  and  especially  with  the  motions  of  both 
the  air  and  rocket  gases  which  create  the  plumes.  It  is  thus  ah  exten¬ 
sion  of  early  work  by  Hill  and  Habert  (1963).  In  the  latter  a  simple 
theory  was  developed  for  the  size  and  shape  of  high  altitude  plumes, 
and  the  dependence  of  these  properties  Oh  altitude,  motor  size,  nozzle 
geometry,  conbustion  pressure  and  forward  speed.  The  general  plume 
was  described  by  a  two -parameter  family  of  ellipsoids;  Using  a  simple 
blast  wave  theory  the  radius  of  curvature  at  the  nose  was  evaluatedj 
while  the  plume  length  was  calculated  using  the  concept  of  plurhe  drag. 
Comparison  between  predictions  and  optical  observations  showed 
excellent  agreement  in  noSe  radii  and  rhaximum  diameters.  In  addition 
it  was  shown  that  it  is  reasonable  to  USe  the  methods  of  continuum 
mechanics  tO  Study  the  flows  in  plumes.  Ah  appropriate  approximation 
to  the  Navieir-StOkes  equations  was  developed  for  the  Structure  Of  the 
mixing  layer  in  the  frOntal  regions,  where  it  is  merged  with  both  the 
external  and  internal  shock  waves.  The  Success  of  the  early  wOrk 
suggested  that  it  wOuld  be  possible  to  develop  methods  to  evaluate  the 
entire  flow  field  in  plumes,  a  goal  Of  great  practical  interest  since  it 
would  permit  calculation  of  air-gas  densities,  chemical  reactiohs  and 
electro -optical  properties. 

This  report  describes  the  application  of  theoretical  principles 
to  the  problem  of  devising  methods  to  calculate  the  Velocities  and 
thermodynamic  properties  of  flow  fields  through  and  arOund  high- 
altitude  plumes.  The  main  effort  is  directed  at  obtaining  a  good  repre¬ 
sentation  of  both  the  jet  flow  from  rocket  motors  and  the  shock-mixing 
layers  in  the  frontal  regions.  Farther  downstream,  and  at  larger 
distances  laterally  from  the  plume  axis,  standard  methods  are  avail¬ 
able.  The  work  may  be  divided  in  three  parts:  the  development  of 
methods  to  calculate  the  flow  fields  throughout  plumes  using  the  principles 
of  inviscid  continuum  mechanics;  the  study  of  the  shock-mixing  layer 
structure  at  the  nose  of  the  plume  using  the  general  Navier -Stokes 
equations  (i.  e.  accounting  for  viscosity,  heat  conduction,  diffusion);  and 
an  estimate  of  the  validity  of  continuum  mechanics  for  describing  high 
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altitude  plumes,  as  inferred  from  a  study  of  shock  wave  formation  at 
high  altitudes. 

It  should  be  noted  that  in  this  wOrk  there  is  the  requirement 
that  the  gross  properties  of  plumes,  such  as  nOse  radius  and  length, 
be  predictable  from  the  detailed  flow  picture,  and  that  they  must  agree 
with  the  simple  rriethods  already  developed  in  the  early  work,  as 
well  as  with  Observation  Of  flights. 
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2.  SUMMARY  ANL»  CONCLUSIONS 


Methods  have  been  developed  to  calculate  flow  fields  in  rocket 
plumes  at  high  altitude.  These,  together  with  the  technique  of 
Witham  (I960)  supply  a  complete  and  practical  way  to  analyze  the  field 
properties  (velocities  and  thermo -dynamic  quantities)  using  the  prin¬ 
ciples  of  inviscid  continuum  mechanics.  The  methods  may  be  des¬ 
cribed  as  follows. 

a.  The  inner  flow  from  the  rocket  nozzle  is  identical  to 
that  which  would  occur  in  a  vacuum.  The  inviscid, 
perfect-fluid  formulation  of  gas  dynamics  is  used  to 
obtain  asymptotic  equations  for  continuous  flows  issu¬ 
ing  from  nozzles  into  a  vacuum.  A  method  of  series 
expansion  is  employed  such  that  the  first  term  repre¬ 
sents  the  final  (limiting)  flow  at  very  large  distances 
from  the  nozzle.  The  general  solution  for  higher  order 
terms  is  demonstrated.  Specific  solutions  for  first 
and  second  order  terms  are  obtained.  The  flow  is 
described  as  a  series  of  perturbations  of  a  flow 
everywhere  radial  and  at  the  limiting  velocity;  the 
angular  distribution  of  mass  flux  along  various  radial 
directions,  according  to  the  theory,  can  be  chosen 
arbitrarily  and  can  therefore  be  used  to  match  the 
initial  conditions  of  the  flow  near  the  nozzle.  The  latter 
conditions  must  be  obtained  independently  by  experi¬ 
ment  or  calculation,  e.  g.  ,  by  the  method  of  character¬ 
istics. 

b.  The  surface  of  the  plume  is  conveniently  described,  in 
the  inviscid  approximation,  by  the  dividing  stream¬ 
line  separating  the  outer  air  flow  from  the  inner  flow 
of  rocket  gases.  Starting  with  the  approximations  of 
hypersonic  flows,  the  inviscid  plume  equation  is  set 
up  which,  when  integrated,  locates  the  surface  of 
dividing  streamlines.  The  equation  is  integrated  for 
representative  cases. 

c.  With  the  dividing  streamlines  as  reference,  the  location 
and  properties  of  the  double  shock  layer  are  expressed 
algebraically. 

d.  Flow  directly  downstream  of  the  rocket  is  returned  to 
atmospheric'  pressure  through  a  mach  disc  which  is 
located  by  the  methods  of  Hill  and  Habert  (1963);  flow 
outside  the  wake  of  rocket  gases  and  behind  the  frontal 
regions  of  the  plume  can  be  analyzed  by  the  technique 
of  Witham.  Thus  the  set  of  methods  (a) -(d)  is  a  complete 
system  for  inviscid  plume  analysis, 
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The  flow  at  the  nose  of  plumes  has  been  analyzed  using  the 
Navier-Stokes  equations  accounting  for  viscosity,  diffusion  and  heat 
conduction.  The  reduced  equations  of  Hill  and  Hubert  (1963)  were  inte¬ 
grated  using  the  method  of  influence  coefficients  to  satisfy  the  split  bound¬ 
ary  conditions.  This  inethod  gave  a  set  of  numerical  results  for  the 
flow  pattern  in  the  merged  layer  around  a  spherical  release;  the  analysis 
represents  the  first  computation  of  its  kind.  It  was  found  however  that 
the  method  was  inherently  awkward  and  slow.  To  remedy  this  a  second 
method  was  developed  for  programming  in  FORTRAN.  The  problem 
is  one  of  integrating  a  set  of  eight  first-order  ordinary  differential 
equations  which  are  restricted  by  a  set  of  split  boundary  conditions. 
Development  of  the  new  method  was  accomplished  through  the  use  of  a 
set  of  asymptotic  (analytic)  solutions  to  the  equations,  an  approach  which 
greatly  simplifies  the  numerical  program. 

The  numerical  results  from  the  first  method  show  the  struc¬ 
ture  of  the  merged  layer  for  one  particular  flight  condition,  and  dem¬ 
onstrate  the  adequacy  Of  the  basic  formulation  of  the  problem  to 
represent  the  variation  o'  properties  in  the  merged  layer. 

The  study  of  shock  wave  formation  at  high  altitude  shows  that 
no  completely  satisfactory  criterion  for  shock  wave  formation  has 
yet  been  established  in  the  field  of  rarefied  gas-dynamics.  A  plausible 
criterion  can  be  developed  however,  expressible  in  terms  of  the  ratio 
of  shock  thickness  to  plume  (or  body)  radius.  Application  of  this 
criterion  to  rocket  plumes  leads  to  the  result  that  typical  ICBMS  and 
IRBMS  will  generate  a  bow  shock  all  the  way  to  burnout  of  the  last 
stage  (SECO). 


4 


iMC  63-80-Rl 


3.  INVISCID  PLUME  FLOW 


3.  1  Getieral 


The  plume  shape  behind  a  rocket  is  the  result  of  interaction 
between  the  atmosphere,  the  forward  motion  of  the  missile,  and  the  jet 
flow  from  the  rocket  motor.  The  inviscid  flow  model  is  that  shown  in 
Fig.  3.  1.  The  coordinate  axes  are  fixed  in  the  missile,  thus  giving  the 
atmospheric  free  stream  a  velocity  relative  to  the  origin.  Behind  the 
rocket  is  a  region  in  which  the  flow  is  identical  to  the  flow  of  a  jet  ex¬ 
hausting  to  vacuum.  This  region  is  bounded  by  an  inner  shock  which 
alters  the  jet  velocity  and  pressure  to  accommodate  them  to  the  on¬ 
coming  free  stream.  The  latter  must  also  accommodate  itself  by  pass¬ 
ing  through  an  outer  shock.  In  between  the  two  shocks  is  a  zone  called 
the  merged  layer  in  which  the  two  gases  flow  together  downstream.  In 
reality  the  gases  mix,  and  the  shocks  extend  well  into  the  layer  because 
of  viscous  forces.  The  properties  of  the  layer  are  complex  and  depend 
on  the  effects  of  viscosi^,  diffusion,  heat  transfer  and  the  properties 
of  real  gases.  A  discussion  and  analysis  of  these  complications  are 
taken  up  in  a  separate  section  (4).  In  the  idealized  inviscid  model  the 
layer  can  be  thought  of  as  flowing  in  two  parallel  sublayers  of  unmixed 
gases,  separated  by  a  contact  surface  as  indicated  by  the  dotted  line  in 
Fig.  3.  1.  Directly  downstream  from  the  rocket,  the  region  of  jet  flow 
is  returned  to  atmospheric  pressure  by  a  combination  of  shock  and  vis¬ 
cous  processes  which  are  not  yet  known.  For  this  model  it  will  suffice 
to  assume  that  the  jet  flow  is  returned  to  ambient  pressure  through  a 
normal  shock  in  the  rear.  This  condition  determines  the  location  of  the 
rear  shock  and  the  flow  behind  it  for  the  inviscid  model. 

In  order  to  analyze  the  inviscid  plume  model  it  is  first  necessary 
to  have  a  knowledge  of  the  jet  flow  from  the  rocket  motor  at  large  dis¬ 
tances  from  the  nozzle  exit.  To  this  end  a  general  analysis  was  made  of 
asymptotic  flows  of  jets  exhausting  a  vacua.  The  results  provide  a 
method  of  representing  the  asymptotic  flow  behind  any  rocket  motor  in 
a  useful  form.  With  this  o.s  a  basis,  the  analysis  of  rocket  plumes  was 
then  carried  out  to  obtain  the  shape  and  position  of  the  contact  surface 
in  the  merged  layer. 

In  Section  3  an  account  is  given  of  the  development  of  general 
formulas  for  the  asymptotic  flows  of  jets.  This  is  followed  by  a  de¬ 
scription  of  the  inviscid  plume  equations  and  their  solution  in  certain 
cases.  These  results  are  then  examined  in  the  light  of  experimental 
data.  Finally  a  discussion  is  given  of  the  next  steps  in  the  analysis  of 
the  rocket  plume  problem. 
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3.  Z  The  Asymptotic  Flows  of  Jets  Exhausting  into  Vacua 

There  exist  simplified  representations  of  flows  from  jets,  such 
as  that  of  Mirels  and  Mullen  (1962)  using  the  concept  of  hypersonic 
similitude,  and  that  of  the  Convair /Astro  group,  using  the  concept  of 
simple  source  flow.  These  methods  aim  at  representing  the  flow  of 
jets  at  large  distances  from  the  nozzle.  There  are  Other  methods  which 
deal  with  the  flow  close  to  the  nozzle,  notably  the  method  of  character¬ 
istics  as  employed  by  Bowyer  (1958).  There  is  a  large  region  of  the  flow, 
however,  which  is  intermediate  between  these  two  zones  of  application, 
and  in  which  either  method  presents  difficulties.  To  carry  the  method 
of  characteristics  downstream  for  thousands  of  nozzle  diameters  is 
impractical.  To  employ  the  source  flow  concept  properly  requires  that 
this  flow  be  matched  to  the  initial  flow  conditions  close  to  the  nozzle, 
but  this  region  is  just  where  the  simple  source  concept  is  not  valid.  The 
same  is  true  of  the  hypersonic  similarity  method.  FOr  these  reasons  a 
new  approach  has  been  taken  in  which  the  rigorous  non-viscOus  equations 
of  motion  are  solved  by  using  an  asymptotic  series  expansion  i.i  terms 
of  reciprocal  powers  of  r,  the  radial  distance  downstream  measured 
from  the  nozzle.  The  first  two  terms  Of  this  expansion  are  equivalent 
to  the  source-flow  concept  of  Mirels  and  Mullen  (1962);  higher  terms  in 
the  series  permit  the  flow  to  be  matched  to  the  flow  close  to  the  nozzle, 
and  give  a  sound  estimate  of  the  errors  involved  in  dropping  higher 
order  terms  from  the  calculations.  A  brief  account  of  the  method  is 
given  in  the  paragraphs  below,  while  a  more  detailed  description  is 
given  in  MITHRAS  report  No.  376A. 

The  basic  equations  for  ihviscid  Compressible  flows  are  de¬ 
veloped  by  Tsien  (1958).  The  natural  coordinate  system  for  this 
investigation  is  spherical,  with  the  origin  located  at  the  rocket  nozzle 
and  the  X3  axis  directed  downstream  along  the  nozzle  axis  (see  Fig.  3.2). 
The  flow  is  rotationally  Symmetric  around  the  x^  axis,  and  has  a  ve¬ 
locity  u  along  the  radius  r  and  a  velocity  v  at  right  angles  to  r.  The 
equations  of  continuity  and  irrotationality  take  the  following  form  in  any 
plane  co  =  const: 


{ 


9v 

ar 


ay 

ae 

(3.  la) 
(3.  lb) 


u  =  radial  velocity 
V  =  angular  velocity 
a  =  local  sonic  velocity 


7 


MC  63-80-Rl 


Coordinate  System  for  Jet  Flow 


The  equations  of  motion  can  be  integrated  independently  of  (3.  1)  to  give 
a  relation  between  the  magnitude  of  the  velocity  and  the  thermodynamic 
properties  of  the  flow; 


[,.  (u/c)^] 

7  =  ratio  of  specific  heats 

c  ~  limiting  velocity  of  flow 


(3.2) 


Equations  (3.  la)  and  (3.  lb)  are  homogeneous  in  dimensional  velocities  as 
well  as  r,  and  may  be  normalized  by  dividing  by  c^  y=!s  where  y  ’'  is  a 
reference  length  which  will  be  taken  to  be  the  throat  radius  of  the  rocket 
nozzle.  All  velocities  can  be  replaced  by  (u/c),  (v/c),  and  r  by  r/y*; 
this  substitution  will  be  assumed  in  what  follows.  The  solution  of  (3.  1) 
by  asymptotic  expansion  is  Carried  out  by  expressing  (u/c)  and  (v/c)  in 
the  form  of  appropriate  series.  These  are: 


u 

c 


c  f 


n  n 


v 

c 


f  £;(9)  r-""*!'-*' 


£„  (9)  =  1 

£'„(e)  =0 

=  [l  -  2n  (7  -  1)] 

(')  =  d/de 

(3.  3) 


The  form  of  v/c  is  not  independent,  but  is  picked  so  that,  together  with 
u/c,  it  satisfies  (3.  lb).  The  series  (3)  are  now  substituted  in  (3.  la)  to 
obtain  the  equations  for  the  f^  which  result  from  setting  the  coefficient 
of  each  power  of  r  equal  to  zero.  The  terms  with  a^  require  the  use  of 
(3.  2).  The  result  of  these  substitutions  is  to  show  that  f ^  (9  )  can  be  an 
arbitrary  function  of  0  while  the  higher  order  f^  are  algebraic  functions 
of  fi  (0)  and  its  derivatives. 


The  thermodynamic  properties  can  be  expressed  in  terms  of  the 
solution  (3.  3)  by  using  the  standard  relationships  of  isentropic  flows. 
Pressure,  density,  and  temperature  in  the  jet  flows  are  given  to  the 
first  order  by; 
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p/p„ 

T/T 

o 


{2  (3 -2t)  (-£,)} 

27 


7 

7-i 


{2  (3-2t)  (-  f,)} 


2 

r 


2  ( 3-27)  (-fj) 
ZfRI) 


1 

7  -  1 


J 


(3.  6) 


To  investigate  the  matching  of  asymptotic  flows  to  near -nozzle 
flows,  it  is  necessary  to  express  some  properties  of  the  flow  to  higher 
Orders  than  the  first.  A  convenient  property  for  this  purpose  is  the 
reciprocal  of  the  square  of  the  Mach  number^  it  contains  in  the 

first  term  the  function  f^  (6)  as  a  linear  factor.  The  near  nozzle  calcu¬ 
lations  of  BoWyer  (1958)  give  (ll^a)  explicitly  and  thus  provide  a  Simple 
way  to  determine  the  form  of  fj  (0).  The  formula  for  a2/u2  =  (I/mZ)  is, 
including  the  second  order  term; 


(7-1)  (3-27)  (-fj) 
“  T(7 -1)- - - 


+  {  3(3-2v)^f^  -2(5  -4v)  •  • 

2r 

(3.  7) 

where 

2(5  -  47)  f^  = 

X  fi  +  S’  fjf'jcote  +  r  f’j  +  T>  fjf'j' 
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1 


and 


=  ~  (3  -  27)^  (27  -  1) 


=  (3  -  27) 

*  4  -  37 

=  (3  -  27) 

To  match  flows  it  is  first  necessary  to  ^'.hobse  a  suitable  form  for  fj(6). 
This  is  most  easily  explained  by  referei.ee  to  the  sketch  of  Fig.  3.3.  The 
flow  emerges  from  the  nozzle  at  the  left  under  known  conditions  of  Mach 
number  4nd  angiUar  distribution.  The  flow  expands  to  some  maximum 
angle  6m*  distributes  itself  in  a  fan-shaped  pattern  of  streamlines 
which  eventually  become  radial.  The  pattern  near  the  nozzle  is  that  cal¬ 
culated  by  Bowyer  (1958);  it  is  straightforward  to  read  off  values  of 
Mach  number  and  compute  l/M^  for  varying  0  at  constant  radii  for 
compai'isoh  with  (3.  7).  For  a  particular  nozzle  a  plot  of  these  >oints 
for  r  S'  60  is  shown  on  Fig.  3.4  (normalized  to  1.0  at  0  ®  O).  From 
(3.  7)  it  is  seen  that  in  the  firfst  approximation  this  curve  is  the  function 
fl  (0).  The  near  nozzle  computations  thus  indicate  the  general  form  of 
fi  (0).  They  suggest  several  possible  curve  fits,  it  will  be  expedient 
to  choose  the  representation 


A 

r 

F 


fl(6) 


X 


X 

(^  “  27  ) 


4  (7  -  1) 


COS0  ~  COS0  M 


1  -  COS0 

effective  maximum  expansion  angle 


(3.8) 


normalizing  constant 


1 


r  2 

(y-  \\ 

5 

7  -  1 

2(1-  cose  j^) 

■V 
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This  curve  fitting  function  is  shown  on  Fig.  3.  4.  The  particular  form 
chosen  here  lends  itself  to  easy  computation  in  the  applications  that 
follow. 


Having  picked  a  suitable  form  for  fj  (6  ),  one  can  now  match 
flows.  For  this  the  nozzle  is  that  of  Figs,  4,  5,  and  7  of  Bbwyer 
(1958),  The  flow  emerges  with  a  semi-divergence  angle  of  15°,  a 
nominal  exit  Mach  number  of  2.  843,  a  nominal  area  ratio  (exit  to 
throat)  of  5:1,  and  a^  of  1.225.  If  the  asymptotic  series  is  a  correct 
representation,  then  the  matching  problem  reduces  to  the  proper 
selection  of  two  parameters:  y'‘\  the  throat  radius;  and  the  nominal 
angle  of  maximum  expansion,  0  \i.  The  first  matching  is  done  along 
the  flow  axis  (0  =  0)  to  obtain  an  effective  value  of  y'*'.  For  this, 
formula  (3.  7)  can  be  reduced  to 


(7-  1)  (  -  V) 
■;2"T7  -  I) 


1  +  I  2(y+  1)  + 


8(7-  1) 

W-  Jy)(\  --c6sl—) 


] 


2r 


(~k) 

^(y-1) 


(3.9) 


Although  this  formula  is  dependent  on  0  j^,  it  is  not  sharply  sensitive  to 
it;  a  value  of  0j^  can  be  picked  based  on  maximum  Prandtl  Meyer 
expansion  angles,  and  then  adjusted  later  without  upsetting  the  first 
calculations. 


The  value  of  y”^  is  at  first  taken  to  be  the  nominal  value,  and 
then  adjusted  slightly  to  obtain  a  curve  fit  at  large  r.  The  result  of 
this  process  is  shown  in  Fig.  3.  5,  in  which  the  solid  line  iS  3.  9  and 
the  dotted  line  is  the  first  term  (the  so-called  linear  solution).  It  is 
seen  that  the  second  order  solution  ten  is  toward  the  linear  as  r  gets 
large,  but  that  convergence  is  not  rapid.  The  exact  calculations  are 
taken  from  Bowyer  (1958)  and  are  shown  as  the  encircled  points.  To 
obtain  a  fit,  a  value  of  y*  was  picked  such  that  v*  equalled  ye.  the  exit 
radius,  divided  by  2.  22;  the  nominal  factor  is  ^5  =2.  24.  This  result 
promotes  confidence  that  the  2nd  order  solution  is  adequate  to  match 
flows  near  the  nozzle,  and  that  the  first  order  solution  is  quite 
accurate  for  r  >  500.  Below  a  value  of  r  =  7,  exact  flow  deviates 
sharply  from  the  asymptotic;  this  point  corresponds  to  the  arrival  at 
the  center  line  of  the  Mach  wave  from  the  lip  of  the  nozzle.  It  is  to  be 
expected  that  such  deviations  would  occur  since  the  asymptotic  flows 
are  assumed  to  issue  from  a  point.  It  should  also  be  noted  that  the 
point  of  origin  for  measuring  r  is  displaced  downstream  from  the 
nozzle.  The  origin  is  found  by  passing  a  line  through  the  nozzle  lip  at 
an  angle  of  »M  and  locating  r  =  0  where  this  line  crosses  the  nozzle  axis. 
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O  •  NUMERICAL  CALCULATIONS  (CONVAIR) 
CD  •  CONVAIR  APPROXIMATION 
0  -  PRESENT  APPROXIMATION,  III  ORO. 


Figure  3.4.  Asymptotic  Jet  Expansion  Theory  Matching  Flows  with 
Near-  Nozzle  Calculations. 
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Figure  3.5.  Asymptotic  Jet  Expansion  Compared  with  Convair's  Machine  Calculations. 


The  next  step  is  to  select  an  effective  value  of  6  This  is 
done  by  exam'ning  the  function  (1/M^)  at  r  =  constant  and  picking  0m 
♦^o  obtain  the  b.st  fit  for  the  variation  of  (I/M^)  with  6.  The  adequacy 
oi  this  procedure  can  then  be  checked  by  comparing  the  flows  at  other 
values  of  r.  In  the  example  here  the  near-nozzle  data  for  r  =  130  were 
used  and  0w  picked  to  be  117.  5°.  The  variation  of  1/M^  with©  at 
r  =  139  is  shownon  Fig.  3.6,  together  with  the  rigorous  near-nozzle 
calculation.  As  a  check,  a  similar  comparison  is  plotted  for  r  =  -12.  It 
appears  that  the  agreement  is  good,  and  that  the  small  changes  observed 
are  adequately  represented  by  the  second  order  theory.  The  linear 
solution,  representing  the  flow  as  r  is  shown  on  Fig  3.  6  as  the 

dotted  line.  It  should  be  noted  that  the  curves  are  normalized  to  equal 
unity  at  0  =  0;  compared  to  the  values  at  the  axis,  then,  the  second 
order  effect  is  to  reduce  (M)^0  =  q  below  its  final  (linear)  value. 
However,  by  referring  to  Fig.  3.  5,  it  is  seen  that  the  second  order 
effect  is  to  increase  the  values  of  (l/M^)  over  the  linear  values  on 
the  axis.  The  two  effects  tend  to  cancel  in  the  off-axis  regions.  Thus 
it  can  be  said  in  the  present  examples  that  the  linear  solution  approximates 
the  exact  solution  to  10  percent  or  better  for  r  =  500,  and  0  ranging 
between  +  100®,  and  that  any  error  will  become  a  successively  smaller 
percent  as  r  increases. 

The  result  of  the  investigation  of  this  section  is  a  set  of  asymptotic 
formulas  which  give  all  flow  properties  of  jets  in  the  regions  of  "inter¬ 
mediate”  flow  and  far  downstream  of  nozzle  exits.  The  method  provides 
estimates  of  errors.  The  validity  and  accuracy  of  the  linearized 
asymptotic  formulas  are  established.  The  pertinent  formulas  are  (3.  3) 
for  velocities,  (3.  8)  for  flow  distribution  and  (3.  6)  for  thermodynamic 
properties. 

In  concluding  this  section,  three  properties  of  jet  flows  are 
developed  w'hich  will  be  needed  in  later  calculation  of  plume  dynamics. 
These  are  the  mass  flow,  total  pressure  force  and  the  total  momentum 
of  jet  flow  across  downstream  boundaries.  The  boundaries  in  question 
are  spherical  (at  constant  r)  and  bounded  by  circles  of  constant  6  ; 
an  example  is  the  right-hand  boundary  of  Fig.  3.  8  .  The  horizontal 
force  at  the  boundary  arising  from  pressure  is  given  by 


F.  =  2 
J 


^•6 

TT  r  \  cos  0-  p^  •  dy 


e 


sin©  COS0  •  p.  d9 


Pj  from  eq.  3.  6 


rn-c.  g  (v.) 


(3.  10) 


g,  h  are  the  functions,  rhj  is  the  total  mass  flow  of  the  jet  and  c- 
its  limiting  velocity. 
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Figure  3.  6  Asymptotic  Jet  Expansion  Comparison  with 

Machine  Calculations  v=  1,225,  0  =  117.5° 

*  m 
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The  form  of  h  and  g  are  easily  determined;  the  main  result 
needed  here  however  is  that  Fj  has  no  zero  order  term,  but  goes  to 
zero  with  r  as  (l/r)2(7~l) .  It  is  understood  that  this  is  the  first  term 
of  an  asymptotic  series,  just  as  are  equations  3.  6. 


The  second  property  of  interest  is  the  integrated  mass  flow  across 
the  same  boundary.  In  the  first  approximation  it  is 


m.(0) 


2  ( 

2irr  J  p.u.  sin©  d0 

0  ^  ^  e  1 

2m.  (-4)  (•^)  J  (-2Xf  sinede 

J  p-  a-=-  ^  1 


r  /cos  0  -  cos0\i\5 

■*ji‘ 'i-r-coss^  )  1 


(3.  11) 


The  mass  flow  between  and  6  through  a  spherical  ring  would  hence 
be  rh.  —  m.{0);  expression  3.  1 1  is  the  first  (zero  order)  term  of  an 
asym]ptotic  Expansion. 


Finally  the  integrated  momentum  across  the  boundary  between  0 
and  0  may  be  computed  as 

r® 

M-(0)  =  J  {dm.(0)/d0|  c.  cos  0  d0  (3.  12) 

J  0  ^  J  ^  ^ 

r  R  ^  ”  cos  0  /  .. 

=  m^Cj  [1  -  g  (0)‘cos0  +  ( - g - £l)(g‘’{e)  _  i)J 

where 

cos  0  -  cos  0,. 

8^®^  ~  1  -  cos  6,  - 

M 

3.  3  The  Inviscid  Plume  Equations 

The  photographic  evidence  of  Rosenberg  (1961)  shows  the  plume 
to  be  roughly  parabolic  at  the  head,  tapering  downstream  to  a  wake 
having  a  maximum  width.  Since  the  forward  speed  of  the  missile  is 
highly  supersonic,  the  flow  field  around  the  plume  is  hypersonic,  at 
least  as  far  back  as  the  region  where  the  hypersonic  parameter  (M^ff*^) 
is  of  the  order  (1).  Here  is  the  free  stream  Mach  number,  and 
a°  is  the  local  slope  of  the  plume  surface  in  radians.  For  a  typical  fcr- 
ward  speed  of  M°  -  7,  •  8®;  this  means  that  the  greater  part  of  the 

head  of  the  plume  is  subject  to  hypersonic  flow  conditions;  this  fact 
will  be  assumed  in  the  initial  analysis.  For  the  purpose  of  setting  up 
force  and  momentum  balances  the  division  between  the  outer  flow  of 
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the  free  stream  and  the  inner  jet  flow  will  be  assumed  to  be  a  surface 
located  approximately  at  the  contact  surface  shown  in  Fig.  3.  1.  The 
first  step  is  to  set  up  a  force -momentum  balance  which  will  determine 
the  location  of  this  surface.  Following  this,  more  detailed  considera¬ 
tions  will  position  the  inner  and  outer  shocks  and  fix  conditions  in  the 
merged  layer.  A  force -momentum  balance  normal  to  the  plume  sur¬ 
face  will  consist  of  three  parts: 

(a)  the  hypersonic  pressure  reaction  as  the  freestreain  enters 
the  merged  layer; 

(b)  the  hypersonic  pressure  reaction  of  the  jet  flow  as  it  enters 
the  same  layer  from  the  inside;  and, 

(c)  the  centrifugal  pressure  of  the  mass  flow  in  the  merged 
layer. 

The  pressures  (a)  and  (b)  will  be  termed  p^^  PML,  J»  respectively. 
The  centrifugal  pressure  (c)  can  be  obtained  from  the  centrifugal  force 
(outward)  per  unit  arc  length  along  the  plume,  designated  AF/AS;  this 
is  the  force  which  would  be  exerted  on  a  strip  of  length  Zirr  sin  6  and  of 
unit  width,  assuming  the  strip  was  straight  so  that  the  pressure  forces 
were  additive.  The  coordinate  system  is  sketched  in  Fig.  3.  7,  and  is 
taken  to  be  consistent  with  that  of  the  previous  section.  The  normal 
force  balance  states  that  (a)  equals  the  sum  of  (b)  +  (c)  at  every  point, 
or 


(3.13) 


(3.14) 


(3.  15) 
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IDEALIZED  PLUME  SURFACE 


Figure  3.  7.  Coordinates  for  Calculation  of  Idealized  Plume 
Surface 
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where 


“  mass  flow  in  the  merged  layer 

V» -T  “  average  velocity  of  the  merged  layer 

R  =  radius  of  curvature  (  +  )  of  the  plume 

surface  in  the  (r,6  )  plane 

Putting  (3.  14)  and  (3.  15)  in  (3.  13,  and  normalizing  by  dividing  by  p, 
(the  chamber  pressure  of  the  rocket  motor),  one  gets  ^ 


(3.16) 


(3.17) 


Equation  (3.  16)  gives  the  inviscid  plume  surface.  It  is  a  second  order 
differential  equation  for  the  spherical  radius  r  as  a  function  of  0.  The 
second  term  on  the  right  accounts  for  the  centrifugal  pressure,  and 
adds  considerable  complexity  to  the  equation.  The  second  derivative 
appears  in  R,  the  radius  of  curvature;  furthermore  the  factor  mML  times 
^ML  is  not  constant,  and  must  be  obtained  from  an  independent  rela¬ 
tion.  In  the  course  of  solving  the  equation  a  simplified  relation  was 
first  used,  obtained  from  (3.  16)  by  dropping  the  centrifugal  term. 

The  results  of  solving  this  equation  did  not  show  a  close  correspondence 
with  the  data.  Also  it  is  possible  to  show  that  the  centrifugal  term 
should  be  retained  on  the  basis  of  an  order  of  magnitude  analysis.  The 
solution  given  here  consequently  will  be  for  the  full  second  order  re¬ 
lation  (3.  16). 

The  solution  of  (3.  16)  requires  expressions  for  (pj/p^.)  Mj  and 
(*^MLVML).  The  first  is  obtained  from  the  jet  flow  functions  of  the 
previous  section  (using  only  the  linear  approximation  as  given  in  eqa. 

3.  6  and  3.  7).  To  obtain  an  expression  for  the  variation  of  (nh 
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one  can  set  up  an  integral  force -momentum  balance  in  the  axial  direction. 
The  control  volume  in  which  the  balance  is  to  be  computed  is  shown  in 
Fig.  3.  8.  The  "ingoing"  quantities  are  free  stream  pressure  and  ■ 
momentum  over  an  area  of  irr^  sin^O  ,  and  the  motor  pressure  and 
momentum.  These  are  balanced  by  the  outgoing  momentum  of  the 
merged  layer,  the  outgoing  momentum  of  the  jet  flow  from  the  right 
hand  boundary  and  the  jet  flow  pressure  over  the  same  boundary. 

The  ingoing  quantities  are 


2  •  2^  ,  o  ,  0..2.  .  2  .  ,  tt2  \ 

wr  sin  elp  +  p  U  )  +  Pex“ex> 

=  ity*^p_,  [2Wq“/p^r)^»in^S  +  C°] 


(3.  18) 


where 


q°  =  ^  p  M  =  free  stream  dynamic  pressure 

=  vacuum  thrust  coefficient 


Vacuum  Thrust 


(p  +  p  ) 

*^ex  '^ex  ex' 


(definition) 


(  )  indicates  conditions  at  rocket  motor  exit. 

'  'ex 

The  approximation  made  in  (3.  18)  is  that  p°/q°  is  small,  i.  e. ,  the  same 
as  that  made  in  the  hypersc  lic  formulas  (3.  14).  The  outgoing  quantities  are 

=  outgoing  momentum  of  the  merged  layer 

M.(0)  =  outgoing  momentum  of  the  jet  flow  across  the 

right  botindary  of  Fig.  3.  12  which  can  be  ob¬ 
tained  from  the  asymptotic  jet  formulas  (See 
Eq.  3.  12). 

Vj^L  "  average  velocity  of  merged  layer. 

”  average  mass  flow  in  merged  layer, 
m.  is  mass  flow  of  rocket  motor  and  is  constant: 

•I  ' 

rhj^L  varies  with  location  along  the  merged  layer. 

F.  =  pressure  force  of  jet  at  right  boundary;  it  can 
^  be  obtained  from  the  asymptotic  formulas  of  (3); 
it  is  to  the  lowest  order  term  (see  Eq.  3.  10): 
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Figure  3.  8.  Control  Volume  Geometry  for  Horizontal  For.e 
Balance 
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ITlC.g(‘Y-)  r  1 


The  expressions  for  Mj(9)^and  Fj  are  the  first  terms  in  the  asymptotic 
series  in  1/  r^’^(7j“l),  Mj(0)  starts  with  a  zero  order  term,  and  Fj 

with  one  of  first  order,  which  may  be  dropped  by  comparison  in  the 
approximation  here.  The  horizontal  force  balance  can  now  be  written 
as 


’’y*^  P<.  [2  Hq‘7 Pj.  sin^e  +  c  p  ] 

and 

Cf  =  Mj(eM>'’'y*^Pc 

The  expression  (rhj^LVj^L)  can  now  be  put  in  terms  of  r,  9  ,  and  a®: 
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p 


where 


cosa 


—  [Zr  sin^9+  \i(y.)g^C 
so  '■  J  \ 


D  COS  H  +  cos  W 


M 


) 


(3.  19) 


r  ='^q°/p_  ‘  (r/y*);  p(v.)  given  below 
^  j 

On  substituting  (3.  19)  in  the  normal  force  equation,  (3.  16),  and  ex¬ 
pressing  R  (the  radius  of  curvature)  in  terms  of  derivatives,  one  gets 
a  second  order  non-linear  equation  for  the  plume  surface.  Using  the 
variables  T,  0  the  equation  can  be  written 


1  dh 
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or  in  functional  notation 
1 


d^T 


d0 


(3.20) 


This  relation  contains,  in  addition  to  the  coordinate  variables  (r,0  ),  only- 
two  explicit  parameters,  (7j»0j^)t  and  these  are  fixed  by  the  gas  constants 
of  the  rocket  motor  exhaust,  and  by  motor  geometry.  The  effects  of 
all  other  parameters,  i.  e. ,  combustion  pressure,  motor  size,  forward 
speed  and  altitude,  are  contained  implicitly  in  T.  From  the  definition 


r  = 


1/2 


{r/y'’')  it  is  seen  that  the  operative  parameter  is  (q°/p^) 


and  that  it  is  simply  a  scaling  factor.  For  a  given  motor  type,  then,  a 
single  computation  will  give  the  plume  shape  for  all  (high)  altitudes,  all 
motor  sizes,  various  chamber  pressures,  and  all  (supersonic)  flight  speeds. 


It  should  be  emphasized  that  (3.  20)  is  based  on  hypersonic  flow 
approximations,  and  therefore  will  not  apply  to  the  rear  portions  of  the 
plume  where  the  slopes  are  of  the  order  of  1/M  or  less.  It  is  also  true  that 
in  the  downstream  region  of  the  plume  the  concept  of  a  layer  contain¬ 
ing  both  exterior  and  interior  shocks  does  not  apply;  the  exterior  shock 
will  become  an  acoustic  wave  propagating  outward  at  the  mach  angle, 
while  the  inner  shock  will  curve  in  and  remain  near  the  axis  (see  Fig  3.  1). 
The  large  region  of  flow  which  is  thus  opened  up  between  them  can  be 
analyzed  using  appropriate  techniques,  which  are  those  of  the  linear  and 
second  order  theories  of  supersonic  flows. 

Solutions  of  (3.  20)  are  presented  in  the  next  section.  In  addition  to 
the  positioning  of  contact  surfaces  which  such  solutions  represent,  it  is 
also  possible  to  derive  simplified  properties  of  the  double  shock  layer 
from  the  inviscid  model.  By  treating  the  layer  as  if  it  were  a  mixture 
of  the  two  gases,  instead  of  being  separated  by  a  contact  surface,  it  is 
possible  to  estimate  the  overall  thickness  of  the  layer  and  its  average 
properties.  To  this  end  an  expression  can  be  written  down  for  mj^jL 
the  mass  flow  in  the  merged  layer,  which  is  the  sum  of  contributions 
from  the  free  stream  air  and  the  jet;  it  is  given  by  (see  Fig  3.  8): 
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The  average  pressure  can  be  expressed  using  the  notation  of  the  normal 
force  balance  (3.  13): 


Pml  Pml.  J 


ML 


which  gives 


^ML  ^  .  2  o 

- =  2  sin  rr 


rh  V 
ML  ML^ 

- Tl - 

Try-  p 


4r  R  sin0 


{3.  25) 


The  thickness  of  the  inviscid  double  shock  layer  may  now  be  estimated. 
Using  6ml  to  indicate  thickness  measured  normal  to  the  plume  surface, 
a  mass  flow  balance  gives 


or 


m 


ML _ 

'ml  ■  ZnrslnH-  P 


Next,  using  an  averaged  equation  of  state 


^zML 


ML 


Pml 

^ML 


the  expression  for  6j^l  rn^y  be  recast  in  terms  of  quantities  already 
obtained: 


-  ^ML 

Zr  sin  0 


(3.26) 


in  which  the  barred  quantities  are  non-dirnensional,  given  by  multiplica¬ 
tion  of  (q°/p^,)  ^'^/y'"'  into  dimensional  quantities.  In  functional  notation 
3.  26  takes  the  form 


This  is  similar  to  the  form  3.  20  except  for  the  appearance  of  a  new 
parameter  (U°/ cj)  the  ratio  of  free  stream  velocity  to  limiting  jet 
velocity.  It  is  interesting  to  note  that  this  parameter  infl>iences  the 
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thickness  of  the  double  shock  layer  but  not,  according  to  3.  20,  its 
position.  Velocity  ratio  (U^/ cj)  is  equivalent,  in  the  hypersonic  approxi¬ 
mation  to  the  root  of  the  ratio  of  stagnation  temperatures;  specifically 


/  c. 
J 


1/2 


An  extension  of  the  analysis  will  give  the  individual  shock  layers, 
i.  e.  ,  the  properties  of  the  air  layer  between  the  contact  surface  and 
the  outer  shock,  and  the  exhaust  gas  layer  on  the  inside.  Assigning 
the  symbol  (  )air  to  the  air  layer  and  (  )jet  to  the  inner  layer,  one  can 
write  four  equations  for  the  two  layer  thicknesses  and  the  two  velocities 
in  the  layers.  They  are 


^ML  "" 

^jet,  ML  ^air,ML 

• 

"'jet,  ML 

ZirrsinG  (*jet ,  ML^^^jet,  ML^  ^^jet,  ML^ 

"'air,  ML 

2Trrsin0{6^.^  j^j^)  ML^ 

"'ml^ml  "" 

^"'jet,  ML^  ^^jet,  ML^  ^  ^"'air,  ML^  ^^air,  ML^ 

The  above  eqxiations  are  (1)  the  total  thickness,  (2)  and  (3)  the  mass 
flow  in  the  air  and  gas  layers,  and  (4)  the  total  momentum.  There 
are  thus  four  eqtiations  for  the  four  unknowns  which  are  soluble  by 
routine  algebra.  The  left  sides  are  all  known  from  the  previous  for¬ 
mulas;  the  densities  on  the  right  are  obtainable  from  the  pressure 
and  enthalpy  by  equations  of  state  of  the  same  form  as  given  for  aver¬ 
aged  quantities  in  the  merged  layer. 

This  completes  the  development  of  methods  for  analysis  of  plume 
flow  in  the  region  of  the  head,  based  on  inviscid  flow  theory.  Numerical 
solutions  and  computations  using  these  methods  are  described  in  the  next 
section. 

3.4  Solution  of  the  Inviscid  Equations 

The  equation  to  be  solved  is  (3.  20),  giving  the  loci  of  plume 
surfaces;  it  is  a  second  order  equation  containing  a  singularity  (  the 
denominator  (c)  of  (3.  20)  passes  through  zero).  The  boundary  con¬ 
ditions  are  two  in  number  and  are  applied  at  the  exit  of  the  rocket  motor. 
The  integration  is  best  carried  out  by  machine  using  standard  procedures; 
that  used  here  is  an  application  of  the  Runge-Kutta  method  expressed 
in  R.  I.  P.  language  on  a  ReComp  III  computer.  Numerical  results  are  best 
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presented  in  graphical  form  showing  loci  of  contact  surfaces,  singu¬ 
larities,  shock  layer  thicknesses,  etc.  Before  passing  to  the  final 
results,  there  are  two  points  to  be  mentioned  which  are  necessary  to 
the  understanding  of  the  integration  and  the  interpretation  of  the  calcu¬ 
lations. 

The  correct  boundary  conditions  must  be  in  accord  with  the 
representation  of  the  jet  by  the  formulas  of  section  3.  2,  The  sketch  of 
Fig.  3.  3  shows  the  situation.  A  necessary  conditic  .  at  the  start  of 
computation  is,  at  6  “  6m»  that  7  sin  0jji  should  equal  (qo/pc)^/^» 
the  curve  start  at  the  lip  of  the  nozzle.  The  second  condition  on  the 
slope  is  arbitrary.  However,  by  definition  a  high  alti  ude  plume  is 
such  that  the  value  of  (qo/Pc)^/^  very  small;  the  details  of  how  the 
contact  iJurface  starts  should  not  be  influential  in  the  final  result. 
Consequently  the  curve  will  be  started  for  0  =  ©m  ^  indefinitely 
small,  a  set  of  conditions  which  give  a  determinate  solution.  The  start 
can  be  approximated  analytically  by  T  =  a  (6m  -  6)”  •  By  substitution 
in  Eq.  (3.  20)  it  can  be  fhown  that  a  and  n  must  have  the  following  values: 


a  = 


3(2  sin^0  ) 
'  m 


(y  sin  0  -  2  ub) 

*  m  ' 


1/2 


where 


b  -  (sin  0  ) /(  1  -  cos  0  ) 

y,  fi,  as  given  in  (3.  20) 
n  =  3 

These  values  can  be  used  to  start  the  machine  computation  at  some 
(small)  value  of  "r  ,  close  to  0  =  ^m* 

The  second  question  is  that  of  the  singularity  in  Eq.  (3.  20) 
when  the  denominator  (c)  is  zero.  Inspection  of  (c)  shows  that  when 
0m  ^  'f/2,  the  value  of  (c)  cannot  be  zero  for  finite  T  sin  0.  However, 
for  ©m  ^  there  is  a  locus  of  points  for  (c)  =  0  which  the  plume- 
curve  must  cross.  By  referring  to  Eq.  (3.  19)  it  is  seen  that  (c)  =  0 
would  make  the  momentum  of  the  merged  layer  equal  to  zero  unless, 
at  the  same  point,  the  slope  c"  of  the  contact  surface  is  vertical. 
Physical  reasoning  thus  requires  that  the  plume  curve  (contact 
surface)  pass  through  the  locus  of  singularities  vertically.  It  can  be 
shov/n  that  the  form  of  Eq.  (3.  20)  is  such  that  if  the  plume  curve 
approaches  the  line  of  singularities  at  an  angle  other  than  rr/Z,  the 
radius  of  curvature  will  rapidly  become  small,  turning  the  curve 
toward  vertical  and  keeping  the  second  derivative  finite;  the  situation 
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has  a  "  static  stability.  ”  In  practice  it  was  found,  in  machine 
calculation,  that  if  the  interval  of  computation  did  not  fall  exactly 
on  the  line  of  singularities,  then  the  computation  proceeded  across 
the  line  vertically  without  special  attention.  The  details  will  be 
illustrated  by  a  specific  example. 

For  specific  computations,  a  single  value  of  Yjet  used: 

Yjet  -  !•  25,  corresponding  to  the  value  for  exhaust  gases  from  a 
lox  -  RP  -  1  motor.  The  values  of  will  be  varied  in  accordance 
with  the  following  scheme:  if  it  is  assumed  that  exhaust  nozzles  have 
parallel  flows  at  the  exit,  then  6m  can  be  taken  as  the  maximum 
Prandtl  -  Meyer  expansion  angle  for  yj  1.  25  and  for  a  Mach  number 
equal  to  the  jet  exit  Mach  number  (Mj  EX)-  Varying  (Mj  EX) 
equivalent  to  varying  0m  (  the  relation  between  Mj  EX  ®m 
shown  in  Fig.  3.  9).  Values  of  Mj  eX  (2, 3,  4,  5)  will  be  used 
corresponding  to  0m  of  (150*  ,  120.5",  99*,  83*,  ).  To  apply 

these  results  to  divergent  nozxles  it  is  only  necessary  to  select  the 
proper  Om  ^^d  interpolate:  the  selection  of  6m  for  rocket  nozzles  is 
described  in  detail  in  section  3.  2. 

The  first  results  are  shown  in  Fig.  3.  10,  which  gives  the  plume 
shape  (contact  surface)  near  the  nose  for  Mj  EX  =  3.  This  shape 
starts  out  at  an  angle  greater  than  ir/2,  becomes  vertical,  then  starts 
back;  the  locus  of  singularities  is  drawn  in  to  illustrate  the  point 
discussed  previously.  The  coordinates  are  non-dimensionalized  by 
iFlo/Pc)^'^/y*>  denoted  by  barred  quantities. 

The  next  figure  (3.  11)  ^hows  a  larger  section  of  the  same  plume 
and  compares  it  with  the  early  calculations  which  neglected  centri¬ 
fugal  pressure  in  the  merged  layer  (equivalent  to  solving  the  first 
order  equation  (b)  s  0,  as  defined  under  Eq.  3.  20).  It  appears  that 
the  simplified  calculations  seriously  tmderestimate  the  size  of  the 
plume  (30  %-  40%  in  terms  of  nose  radius);  it  also  can  be  shown 
that  the  behavior  of  the  contact  surface  as  0  decreases  is  basically 
different  with  the  centrifugal  term.  It  is  concluded  that  the  centri¬ 
fugal  effects  must  be  accounted  for  to  obtain  realistic  results. 

Fig. 3.  12  shows  a  family  of  plume  contact  surfaces  for  yj  =  1.  25. 
A  one  parameter  family  is  si^icient  to  represent  all  possible  rocket 
plumes,  at  least  in  the  forward  regions.  The  parameter  Om  ts 
tabulated  on  Fig. 3.  12,  as  well  as  the  corresponding  value  of  the  exit 
Mach  number  from  a  non-divergent  nozzle.  The  contact  surfaces  of 
Fig.  3.  12  may  be  used  as  the  bases  for  investigation  of  the  thickness 
of  the  inviscid  shock  layers,  the  effects  of  heat  conduction  and 
viscosity,  and  the  continuation  of  the  plume  flow  downstream  of  the 
head  wave. 
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Figure  3.9-  Prandtl -Meyer  Turning  Angles 
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Figure  3.  10.  Plume  Detail  Near  Origin 
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Figure  3.  12.  Plume  Contact  Surfaces  y.  =  1.  25 
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3.  5  Comparison  With  Flight  Data 

The  technique  ofRosenberg  et  al  (.1961)  provides  scaled  photo¬ 
graphs  of  the  chemiluminescence  given  off  by  rocket  plumes  at  night. 
Assuming  that  the  outer  edge  of  the  luminescence  defines  the  outer  edge 
of  the  gaseous  plume,  the  nose  radii  as  determined  by  Rosenberg 
may  be  compared  with  those  obtained  from  the  calculations.  It  is  assumed 
that  the  rocket  motor  can  be  represented  by  the  parameter  =  4 

and  Yj  =  1.  25,  the  nominal  values  for  an  Atlas  sustainer  engine. 

The  effective  nose  radii  of  the  calculated  plumes  are  determined 
in  a  manner  similar  to  that  used  in  obtaining  nose  radii  from  photographs. 
The  nose  is  assumed  to  be  roughly  equivalent  to  a  parabola;  the  parabola 
constant,  which  is  twice  the  nose  radius,  is  evaluated  by  replotUng  the 
curves  of  Fig.  3.  12  in  terms  o^  y^  vs  x.  The  slopes  are  then  the  para¬ 
bola  constants,  as  shown  in  Fig.  3.  13.  Using  the  value  of  R  for 
Mjex  =  4,  the  nose  radius  is  given  theoretically  by 


o 

Y 


Yj 

M. 

jex 


1.4 
1.  25 
4 


This  relationship  is  plotted  on  Fig.  3.  14  together  with  the  data  available 
from  flight.  Included  among  these  data  are  several  flights  and  two  dif¬ 
ferent  missiles.  Also  shown  is  the  curve  from  the  blast  wave  theory 
of  Hill  and  Haber t  (1963).  It  is  seen  that  both  theories  approximate 
the  data  adequately,  particularly  in  view  of  the  uncertainties  involved 
in  the  data  reduction. 


3.6  Discussion 


It  appears  that  the  methods  of  flow  field  analysis  developed  in 
Section  3  are  capable  of  providing  plume  shapes  in  agreement  with 
the  known  data  and  with  the  blast  wave  theory  of  Hill  and  Habert.  The 
flow  field  technique  provides,  in  addition,  formulas  for  the  distribution 
of  thermodynamic  properties  in  the  air  flow,  the  jet  flow,  and  the 
merged  layer  (as  far  as  inviscid  mechanics  is  applicable).  The  de¬ 
velopment  here  has  concenti*ated  on  the  frontal  regions  of  the  plume, 
since  this  is  the  initial  disturbance  in  the  creation  of  the  flow  pattern. 
There  are  standard  methods  avHlabie  for  continuing  the  calculation 
outward  and  rearward  from  the  frontal  regions,  as  for  example  those 
of  Witham  (1950).  The  entire  flow  near  the  axis  (trajectory)  is  thus 
subject  to  analysis  by  these  methods. 


35 


MC  63-80-RI 


Figure  3.  13.  Parobola  Constants 
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There  will  be  modificatioD.ij  of  such  calculations  which  account 
for  the  influence  and  lixTiits  of  some  of  the  assumptions  made  here.  The 
estimation  of  viscous  effects  ^  lust  be  based  on  studies  of  the  Navier- 
Stokes  equations  as  applied  to  this  problem;  Section  4  is  a  detailed 
investigation  into  dissipative  effects.  Another  modifying  condition  is 
the  stratification  of  the  thermodynamic,  chemical  and  electric  proper¬ 
ties  of  the  atmosphere  and  their  effects  on  plume  flows.  In  addition 
the  influence  of  the  earth's  magnetic  field  on  the  electrified  particles 
in  the  flow  become.s  relatively  larger  with  increasing  height.  Finally, 
proper  application  of  the  theory  requires  some  estimation  of  the  limit 
of  validity  of  continuum  mechanics  which  underlie,  the  present  approach. 
This  matter  is  taken  up  in  Section  5. 


4.  MIXING  LAYER  STRUCTURE  AT  THE  NOSE  OF  THE  PLUME 


4.  1  General 


One  important  source  of  observables  associated  with  the  high- 
altitude  rocket  plume  is  the  mixing  of  the  exhaust  gases  with  the  ambient. 
This  mixing  occurs  in  a  layer  whose  mean  surface  coincides  approxi¬ 
mately  with  the  contact  surface  of  the  inviscid  analysis.  In  this  layer 
the  effects  of  diffusion,  viscosity  and  heat  transfer  determine  the 
character  and  properties  of  the  flow.  The  thickness  of  the  mixing 
layer  is  fixed  largely  by  the  Reynolds  number  and  a  typical  plume 
dimension  such  as  the  radius  of  curvature  at  the  nose.  Above  100  Km 
the  Reynolds  number  is  usually  less  than  100  and  the  mixing  layer 
fills  the  entire  region  between  the  inner  and  outer  shocks  of  Fig.  3.  1., 
at  least  for  an  appreciable  distance  back  from  the  nose. 

Galcula  'ons  of  mixing,  layers  in  this  flight  regime  near  the 
limit  of  contints  ri  flow  have  not  previously  been  published;  the  follow¬ 
ing  analysis  represents  the  first  computation  of  its  kind.  It  is  related 
to  several  o'  "merged-layer”  analyses  of  the  flow  about  solid  spheres 
(Levinsky  a:.<  .oshihara,  19ol;  Kb.o,  1963),  but -there  are  several  essen¬ 
tial  differences  due  tc  the  change  in  the  inner  boundary  conditions. 

In  the  first  analysi''  is  necessary  to  restrict  the  study  to  the  region 
of  the  nose  because  u*  mathematical  complexities. 

Sections  4.  2  and  4.  3  present  the  equations  and  approximations 
Avhich  are  employed  to  describe  the  flow  in  the  merged  layer,  and 
explain  a  rhethod  which  has  been  developed  for  solution  by  numerical 
integration  which  exploits  a  set  of  asymptotic  (analytical)  solutions. 

This  method  has  been  selected  on  the  basis  of  experience  with  several 
methods,  and  is  felt  to  provide  the  simplest  and  most  easily  applied 
program  for  obtaining  solutions.  In  the  process  of  developing  numeri¬ 
cal  methods,  several  sample  solutions  have  been  obtained  which  demon¬ 
strate  the  adequacy  of  the  equations  to  describe  the  two  thick  shocks 
merged  with  the  mixing  zone.  In  Section  4.  4  sample  results  are  pre¬ 
sented  which  were  obtained  using  an  earlier  method. 
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4.  2  FORMULATION  OF  THE  PROBLEM 


The  spherical  system  of  reference  chosen  is  depicted  in  Fig.  4.  1. 
The  velocity  components  are  shown  in  the  positive  directions.  The  flow 
is  axisymetric.  The  following  simplifying  assumptions  concerning  the 
fluids  have  been  adopted: 

(aj  The  two  fluids  are  two  perfect  chemically  non-reacting 
gases,  with  constant  specific  heats 

(b)  They  diffuse  mutually  following  Fick's  law,  thus; 

p'v-Vc  -  JtfvP'^c  ~  =  0 

(pI  The  Lewis  number  is  unity,  and  thus  the  expression  for 
the  heat,  transfer  vector  is  not  changed  by  the  mutual 
diffusion  of  the  two  gases. 

(d)  The  Prandtl  number  is  ^ 

(e)  The  viscosity  coefficient  is  proportional  to  the  square  root 
of  the  temperature 

4.  2.  1  Governing  Equations 

The  varuiblcs  describing  the  flow  field  are  non-dimensionalized 
by  reference  to  the  upstream  conditions  according  to  the  following  manner: 


p 

.  V 

u 

„  H 

■  ^ 

V  -  _  u 

00 

00 

^  =  TT 
7V00 

T 

=  X 

Po  '  —7:2 

P2  •  “TIZ 

(4.  1) 

Xi 

8 

< 

8 

^00  00 

where  the  symbols  with  the  subscript  denote  upstream  quantities,  and 
the  overlined  symbols  denote  physical  variables. 

The  set  of  governing  equations  are: 

(a)  mass  conservation  equation 

(b,  c)  radial  and  tangential  Navier  Stokes  equations 

(d)  diffusion  equation 
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INSIDE  BOUNDARY 


OUTSIDE  BOUNDARY 


Figure  4.  1.  Reference  System 
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(e)  energy  equation 

(f)  equation  of  state 


(g)  viscosity  law 

Explicitly  in  spherical  coordinates  they  are: 
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-  'r  -  I 

-  (v  +  u  cot5)j 

^  ^  <7)  7  I7  ^  7  ^ 


(v  +  u  cot©)*  +  r^ 


+  ^  (f>  ^  * 


_1  9u  Svl 

7  F7 


.  1  4  r  8  f.  8hl  .  9  ,M-  3h  ,  _^  8h  ^  cot8  8h‘ 


3}^ 


where 


«  9v  .  1  9u  ,  2v  ,  u  ,rt 

^  ■^73F'^~‘^7 


H  =  h  +  v^  +  u^ 


(4.  2  Coht) 


f  [  P  =  rph 
where 


,  R  .  +  (R  -  R  .  )  c 
p  _  1  air  gas  air 

+  (C 


7  zr 


-C  ■■); 


^air  ^gas  ^alr 
_  T  tl  _  /  h  .-y 

4.  2.  2  Separation  of  Independent  Variables 

The  present  method  of  separating  the  independent  variables  has 
been  used  by  several  investigators  (Hays  1959»  Levinsky  1961,  Kao  1963): 

V  ~  V  cos  9 

o 

u  ~  sin0 

P  =  P« 


u  =  p  cos  6 
o 
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h  cos  5 
o 


h 

H  =  H  cos“0 

o 

p  =  p  cos^d  +  p  sin^d  (4.  3) 

O  to 

The  quantities  with  a  subscriptare  functions  of  r  only, 

4.  2.  3  Derivation  of  the  Reduced  Equations 

In  the  system  of  governing  equations  the  terms  of  order  smaller 
than  A/Ror  p/Re  are  neglected  (A  is  the  thickness  of  the  layer). 

The  expressions  assumed  for  the  separation  of  independent  var¬ 
iables  are  substituted  into  the  full  system  and  yield  the  following  set  of 
reduced  equations,  the  only  independent  variable  being  r. 


(Pv)' 

T 

(u  +  v)  =  0 

Continuity 

pvv' 

+  Po 

=  ^  <T 

Rad.  Mom(l) 

Pz  = 

^-(u  +  v) 
r 

Rad.  Mom(2) 

i 

1 

pvu' 

+ 

r 

(u  +  v)  +  1  (p 

,  1 

(pu’)' 

Tan.  Mom 

1 

pvc' 

1 

Re 

7  (pc')' 

Diffusion 

! 

I 

} 

pvH' 

1 

7  (NtH-)' 

Energy 

1 

I 

i 

Po  = 

U.  — 

rph 

h 

1 

7 

1st  State 

Viscosity 

1 

2G 

T 

00 


00 


+  <i-c)  -r 

Y  ^  M 
00 


H  =  h  +  v2 


(4.4) 
Tot.  Enthalpy 


As  a  matter  of  convenience  the  variables  appearing  in(4v3)are  denoted 
without  indexes  when  the  meaning  of  the  symbol  can  be  determined 
without  ambiguity.  This  convention  is  carried  on  henceforth. 
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-yife,  -J 


The  boundary  conditions  in  the  inner  edge  are  given  by  the  inner 
source  flow.  The  inviscid  system  of  equations  is: 

(pv)'  +  (u  +  v'  =0 

pvv'  +  =  0 

_  pu(u  -f  v) 

^2  r  (4.6) 

Pq  =  r  P  (H  -  v^) 

I  P  u  (u  +  v)  ,  2  ,  A  ft 

pvu'  +  a — i- - L  +  -  (p^  -  p^)  -  0 


H'  =  0 

c.'  '  “  0 

This  yields  rigorously: 

H  = 

c  =  Constant 


Constant 


(4.  7) 
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The  expression  of  the  third  constant  deserves  some  attention. 
The  entropy  of  a  perfect  gas  is  given  by  the  differential  expression: 


dS 


per 


■  1 
y  ^  1 


dv 


where  S  and  v  designate  entropy  and  specific  volume  respectively.  dS 
is  independent  of  the  mode  of  non  dimensionalization,  and  hence: 


dS 


R 


per 


[y-n- 


dh  ,  dv 

T  ^  — 


or 


Sq  being  an  arbitrary  constant  which  will  be  taken  as  zero.  The  appli¬ 
cation  of  this  definition  at  the  inside  boundary  yields  a  physical  meaning 
to  the  third  constant; 


where  Sq  is  the  enthalpy  of  the  gas.  The  two  constants  Hq  and  Sq  are 
provided  by  the  characteristic  features  of  the  source.  The  specific 
application  of  these  conditions  to  the  numerical  solutions  will  be  made 
in  section4.  3.  4.  It  should  be  noted  that  the  last  condition  is  equivalent 
to  saying  that  the  flow  becomes  isentropic  along  the  axis  of  flow. 


P 

— mz 

(H-v^) 
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4.  3  SOLUTION  OF  THE  EQUATIONS 


4.  3.  I  Reduction  to  First  Order  Equations 
By  introducing  the  new  variables 

0  =  jiH';  W  =  ^iv';  x  =  u  =  iiu' 


the  system  of  reduced  equa.tions  of  the  last  section  may  be  recast  as 
the  following  system  of  ten  first  order  equations: 


vp' 

-p. 

2(u  +  v)  ^  W  1 

1 

r  p 

v' 

= 

w/p 

u' 

(A)/p 

C 

= 

x/p 

Pz 

= 

pu(u  +  v)/r 

X  ' 

= 

(3/4)Re  (pvx)/p 

w 

(3/4)Re  (p^  +  pvv') 

w' 

= 

Re" 

|p  [(vci))/p  +  u(u  +  v)  /rj  + 

and 

V. 

H' 

o 

= 

= 

(3/4)Re(p  v^)/p 

(4.8) 


The  ten  unknown  variables  appear  as  first  derivatives  on  the  left  side  of 
the  equations.  The  quantities  p  and  po  on  the  right  side  are  not  to  be  con¬ 
sidered  as  additional  luiknowns,  but  as  functions  of  the  ten  primary  un¬ 
knowns;  they  are  given  algebraically  in  terms  of  the  latter  by 


Po 


where 


h 


o 


r 


rpho 

h 

O 

IC  T~ 

p  00  _ 

-  G  +  (1-  C)  - -^-r 


(4.9) 


,  C  R  +  (1-C)  R  . 

1  gas  '  '  air 

2  CTT — +  (i-cj  r  " 

Pgas  Pair 
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In  the  .gC.^eral  case  this  system  of  equations  requires  specification  of  ten 
constants  of  integration.  If  these  can  be  determined  at  one  point  through 
a  knowledge  of  the  values  of  the  dependent  variables  at  that  point,  then  a 
solution  of  the  equations  is ^  in  principle,  determined.  The  problem 
treated  here  however  has  some  conditions  given  at  the  outer  botindary 
and  the  reat  at  the  inner  boxmdary.  The  handling  of  this'  two-ppint  boun¬ 
dary  value  problem  requires  special  attention,  and  is  discussed  in  the 
next  section. 

The  ten  eqviations  may  be  reduced  to  eight  by  extracting  two  inte¬ 
grals  from  the  system.  The  equations  for  X '  and  ^ '  may  be  combined  to 
give 

=  X'/  X 

so 

^  =  kj  X 

iJsing  the  equations  for  C  and  H',  this  last  becomes 

H'  *  k,  C» 
o  1 
or 

=  kjG  +  k2  (4.  10) 

This  relation  may  be  considered  as  another  algebraic  equiyale-'cc  of  the 
type  (4.  9):  with  it  Hg  may  be  eliminated  as  an  unknown.  The  equation  for 
may  be  dropped,  as  well  as  equation  for  ^  '  since;  ^  appears  only 
in  the  equations  for  H'  and  ^  This  leaves  eight  equations  (the  first 
eight  of  the  set  (4.  8),  and  eight  unknowns  (appearing  as  first  derivatives 
in  the  left  Column),  and  eight  arbitrary  constants  of  integration.  The  two 
constants  ki  and  k2  of  (4.  10)  are  arbitrary  and  complete  the  total  of  ten. 
It  should  be  noted  that  the  linear  dependence  of  total  energy  on  concen¬ 
tration  is  a  special  result  of  picking  the  Prandtl  number  to  be  3/4. 


The  values  of  k^  and  k2  are  determined  from  two  of  the  boundary 
conditions  given  in  the  previous  section.  One  condition  is  at  the  outer 
boundary,  the  other  at  the  inner.  Specifically 


where 


Hr  - 

vj 


1 


(4.  11) 


ratio  of  specific  heats  in  the  ambient 

total  enthalpy  of  irnier  flow;  this  is  the  constant 
of  the  first  equation  in  (4.  7). 
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The  integration  of  the  eight  equations  can  now  proceed  providing 
a  method  is  adopted  which  overcomes  the  problem  of  two -point,  boundary 
values.  A  second  problem  will  be  encountered  at  the  stagnation  point 
(where  v  becomes  zero)  owing  to  the  riature  of  the  first  equatir^n  of  the 
set  (4.  8).  It  can  be  seen  that  the  value  of  p'  will  become  either  infinite 
or  indeterminate. 


The  two  point  boundary  problem  will  be  dealt  with  by  using  the 
method  of  asymptotic  solutions  at  the  outer  boundary,  and  is  taken  up 
in  the  next  section.  The  solution  near  stagnation  will  be  obtained  by 
series  expansion*  and  is  described  in  section  4.  3.  3. 

4.  3.  2  Asymptotic  Solutions  at  the  Outer  Boundary 

The  method  of  asymptotic  solutions  is  to  linearize  the  equations 
at  the  outer  boundary,  where  the  unknowns  may  be  expressed  as  small 
perturbations  from  the  limiting  bounda.ry  values.  The  resulting  equa¬ 
tions  can  then  be  solved  by  standard  methods.  The  solutions  will  include 
ten.  arbitrary  constants  corresponding  to  the  ten  boundary  conditions 
of  the  problem.  As  will  be  shown,  seven  can  be  imposed  in  advance  on 
the  basis  of  boundary  values,  while  another  can  be  fixed  in  advance  by 
relating  it  to  the  Reynolds  number  selected  for  the  computation.  The 
remaining  two  can  only  be  picked  by  relating  them  to  the  inner  flow 
through  numerical  integration  across  the  merged  layer.  How  this  may 
be  done  is  detailed  in  this  and  the  next  section. 


The  asymptotic  forms  may  be  written  as  follows,  using  A  as  a 
modifier  to  indicate  small  quantities: 


p  -  1 +  Ap 

P2  =  ^2 

V  =  -1  +  Av 

X  =  A  X 

u  =  1  +  Au 

W  =  AW 

O 

II 

u  =  A « 

(4. 12) 

These  forms  are  substituted  in  the  set  of  eight  equations  (4.  8); 

eight 

linear  equations  are 

obtained  by  rejecting  terms  of  second  and  higher 

orders.  Theresult 

is 

Apf 

rz(iSu+A.')  ^ 

(4.13) 

Av' 

AW;  AC'  =  AX;  AX'=  -  (l)Re-A  X 

An'  = 

Aw;  Ap^  =  ;  AW'  =  (|-)  Re(Ap'  AW) 

Ad  = 

Re  (-  Ae.  +  ^  -  /Ip  ) 

r  r 

b 

► 
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The  linearized  algebraic  relations  cor^  esponding  to  (4.  9)  are 


A-^C  4 


A  =  < 


/di’\ 

y  c  =  0 


A 

2 


'ir.\  .  l^Bgas 

/ C  =  0  2  c 


air 


(4.  14) 


kj  from  eq.  4.  1 1 


+ 


1  4  Ap  ,  Ap  not  needed  . 


Now  the  equations  for  AX'  and  AC'  can  be  solved  independently;  they 
yield 

AX  =  H3/4)Re  k^e  -<3 /4)Re  r 

AC  =  j^AXdr  =  j^^e“(3/4)  Re  r  ^  (4.  k5) 

where  in.  -  0 
4  • 

Another  integration  is  possible  by  combining  equations  for  p',  v*  and  p^'; 


4Ap' 

^2 


22^^  4  Av’ 

Ap  -  Av 

2  ^  5 


where  k^  =  0 


(4.16) 


These  results  may  be  used  to  reduce  the  set  (4.  13)  to  three  equations  for 
three  unknowns.  It  is  convenient  at  the  same  time  to  redefine  the  in¬ 
dependent  variable  to  be  r  ~  r  Re  and  the  symbol  for  differentiation  to 
be  d/d  r  =  (  )V  s  ()•  /Re.  The  three  equations  for  Au,  Av  and  APare 


(.a) 

(b) 

(c) 


A  p^  -  Av^ 


Av 


W. 


Att^  4 


2(Au  4-  Av) 

T 

I 


1 


=  ~  <Au  4  Ap  -  2Ap 


(4.  17) 
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It  should  be  noted  that  Re  does  not  appear  explicitly  in  this  set.  These 
relations,  although  linear  and  comparatively  simple,  cannot  be  solved 
exactly  in  closed  for.m.  Useful  solutions  are  obtained  by  noting  that  in 
the  regions  of  interest  we  have  t  >10,  and  that  what  is  needed  is  the  be-r 
havior  of  u,  v  arid  p  as  r-*  This  suggests  a  series  expansion  of  the 
solution  in  reciprocal  powers  of  t,  having  the  following  form: 


S 

n“o 


n 


In  this  expression,  the  [^In  are  functions  of  r,  and  will  turn  out  to  be 
solutions  of  linear  equations  with  constant  coefficients.  The  power  p  is 
included  to  account  for  the  fact  that  some  solutions  are  not  expressible 
as  integral  powers  of  (l/r),  and  require  a  fractional  power;  p  is  taken  to 
lie  between  0  and  I.  Its  specific  value  i-s  determined  when  is  substi- 
tutedin:(4.  17),  together  with..similar  expressions  for  ^  and  Av.  Because 
of  the  large  values,  of  t  it  will  be  sufficient  to  carry  the  determination  of 
the  n  only  so  far  as  is  necessary  to  bring  in  the  arbitrary  constants 
of  integration;  thus  will  turn  out  to  be  only  the  zero  order  terms  for  Ap 
and  Av,  and  the  zero  and  first  order  for  Au.  Inspection  of  the  sot  (4.  17) 
shows  th3.t  there  will  be  five  arbitrary  constants  of  integration,  pne  is 
deternraned  at  once  ^by  single  hvtegration  of  (b);  calling  it  k  o>  it  is  seen 
at  once  that  =  0  is  required  by  the  boundary  conditions.  Of  the  re¬ 
maining  four,  two  will  hot  be  determined  by  the  outer  boundary  condi¬ 
tions;  the  significant  solutions  of  (4.  17)  will  contain  these  two  arbitrary 
constants.  The  equations  for  the  zero  order  terms  result  from  direct 
substitution,  and  are  (for  any  value  Of  p) 


E 

o  <1 

1 

N.O  =  "■ 

Mo  =  1  {Mlo  -  M„: 

(<[^o]o 

and  AC  from  4-14  and  4-15) 

These  expressions  have  constant  coefficients  and.  their  solutions  are 
easily  obtained.  They  are 


M]o  " 

I^'^]  0  ‘^7 

Mo  = 

3b2 

k  e  “T  ^  —  - 

^8 

Mo  = 

k,„ 

m2 --  1 

00 

B 

_  - - 

(4.  18) 


where 
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Each  of  these  terms  has  yet  to  be  multiplied  by  its  respective 
To  determine  the  jp*s  it  is  necessary  to  examine  the  higher  order 
terms;  it  is  found  that  the  solution  cannot  be  continued  unless,  for  the 
.k,9  term,  p  s  1;  and  for  the  kg,  term,  p  «  2/  (M^  -  1).  For  the  coi^stant 
terms,  k?  and  k  _  there  are  no  values  of  p  which  will  permit  the  solution 
to  be  continued  in  a  convergent  series,  and  hence  they  are  zero.  This 
completes  the  information  needed  for  the  asymptotic  solutions,  Avhich  can 
now  be  written  out  explicitly; 


Av 


*=3^  ^  -(3/4).t  ^  1 

?ri  rP  '  ■■■ 


M  >  'v/3 


There  are  three  constants  (  k^,  k^  and  k.)  of  integration  which  are  still 
undetermined.  The  seven.already  determined  are  essentially  fixed  by 
the  six  independent  conditions  at  the  outer  boundary  and  one  (total  enthalpy) 
of  the  inner  flow.  Thus  k^,  kg»  and  k„  are  related  to  the  as  yet  un> 
specified  character  of  the  inner  flow.  \i>ne  condition  can  be  applied  as 
follows:  select  a  Reynolds  niunber  based  on  the  distance  from  the  origin 
to  the  outer  boundary  as  determined,  eay,  by  the  value  of  r  at  that  point 
where  the  density  is  1%  above  free  stream  value;  then  setting  Re  =  r 
and  Ap  =  .  01  in  the  asymptotic  formula,  a  condition  is  obtained  on 
and  1^.  This  condition  essentially  fixes  Uie  scale  of  the  flow;  with 
it  one  of  the  Constants  say  kg,  can  be  eliminated.  The  two  remaining 
conditions  are  obtainable  only  by  picking  trial  values  of  the  k.3  and  k9, 
integrating  the  equations  numerically  through  the  merged  layer  and  thus 
relating  the  constants  in  the  asymptotic  solution  to  the  properties  of  the 
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inner  flow.  To  match  a  given  inner  flow  it  will  be  necessary  to  choose 
two  constants  at  the  outer  edge  which  give  two  specified  properties  of 
the  inner  flow  (in  addition  to  the  stagnation  enthalpy,  which  i.s  already 
determined).  This  process  is  essentially  the  inversion  of  a  rank  two 
matrix  and  presents  no  real  problem  once  a  series  of  integrations  has 
been  carried  out.  A  discussion  of  the  integration  through  the  stagnation 
point  and  the  inner  flow  properties  follows. 

4.  3.  3  The  Solution  Near  Stagnation 


As  noted  in  the  previous  secticn,  the  expression  for  p'  in  the  set 
4.  8  becomes  sin/^ular  when  v  =  0.  In  numerical  integration  by  machine 
the  consequence  is  an  instability  of  computation  near  stagnation.  A 
method  of  extrapolating  across  the  unstable  band  may  be  developed  based 
oh  a  series  expansion  about  this  point.  It  is  first  useful  to  select  a 
point  in  the  computation  which  is  near  the  point  of  v  =  0  but  before  insta¬ 
bility  appears,  and  this  requires  a  criteria.  From  the  first  equationof  (4.  8) 
it  is  seen  that  both  v  and  the  bracket  on  the  right  (Should  pass  through  zero 
linearly  with  r,  and  at  the  same  point.  Using  a  scheme  of  linear  extra¬ 
polation,  the  point  where  v  =  0  is  estimated  by 


where 

so 


0 

0 

r 

n 


V  +  v'lf? 

wn 


V  + 


14 

+  rj 


pv 


r^  =  point  at  which  v  =  0 


The  value  of  the  right  hand  bracket  may  be  estimated  similarly  for  the  same 
point; 


The  derivatives  on  the  left  are  all  available  from  the  computations.  The 
value  of  €  should  go  tc  ^ero  faster  than  v  does.  In  practice  however,  it 
may  not.  If  this  is  owing  to  computational  instability  then  the  computation 
should  be  rejected  after  c  has  passed  through  a  small  minimum,  thus 
establishing  the  band-width  across  which  the  series  expansion  is  to  be 
used.  If  on  the  other  hand  c  does  not  pass  through  zero  or  a  small  mini¬ 
mum  for  a  given  pair  of  values  of  k.3  and  kq,  then  these  must  be  adjusted 
to  give  pairs  of  values  for  which  e  does  become  small,  this  process 
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being  essentially  a  restriction  on  the  initial  constants.  Having  arrived 
at  values  of  r)  and  r^^  a  series  expansion  may  be  employed.  Let: 

V  =  (v)^  +  (v)jn  +  .  . . 

u  =  (u)^  +  +  . .  etc. 

All  quantities  in  (  )  are  constant  coefficients. 

Substituting  these  in  the  set  (4.  8)  letting  (v)^  ~  0,  the  coefficients 
of  the  first  order  (first  derivatives)  may  be  determined  in  terms  of  the 
zero  order  coefficients.  The  selection  of  both  v  and  the  right  side  of  the 
first  equation  to  be  0  at  Tq  is  equivalent  to  two  conditions;  there  are 
eight  remaining  zero  order  coefficients  which  can  be  independently 
picked,  giving  a  total  of  ten  conditions  as  before.  The  first  order  coef¬ 
ficients  thus  determined  will  be  substantially  the  same  as  the  first  de¬ 
rivatives  as  computed  from  (4-  8)  at  the  edge  of  the  instability  band,  with 
the  possible  exception  of  the  expression  for  (p)i,  which  is  given  by  the 
expression 
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This  formula  for  (p)^  is  non-singular  and  is  given  in  terms  of  the  zero 
order  coefficients  at  me  stagnation  point.  If,  however  the  values  of  u, 
p  etc.  ,  are  used  from  machine  results  at  the  outer  edge  of  the  instability 
band,  then  (for  small  rj)  the  error  committed  will  be  of  higher  order  than 
the  approximations  given  here.  Thus  a  stable  value  of  dp/ dr  is  obtained 
whicii  is  consistent  with  the  conditions  at  stagnation.  Of  course  the  values 

of  the  quantities  (u,  p - ).  and  their  first  derivatives  may  now  be  used  for 

a  linear  extrapolation  across  the  instability  band  to  re-establish  the  ma¬ 
chine  computations  with  correct  values  of  the  dependent  variables.  The 
integration  may  continue  to  the  inner  flow  and  be  matched  as  described 
in  the  next  section.  For  completeness  the  remaining  relations  are  given 
here  for  the  series  solutions  at  stagnation: 
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4.  3.  4  Matching  Solutions  With  the  Inner  Flow 

The  discussion  of  inner  boundary  conditions  in  section  4.  2.  4  shows 
that  when  inviscid  flow  is  reached,  the  equations  will  automatically 
require  the  zero  order  enthalpy  Hq  to  be  a  constant  (Hq);  sirnilarly 
the  value  of  C  will  becorne  constant.  The  value  of  Hq.  can  be  fixed 
a  priori  by  proper  choice  of  kj,  as  given  by  formula  (4-11).  Aceprding 
to  |4.  7i),  then,  there  remain  two  disposable  properties  of  the  inner  flow: 
the  value  of  C  and  of  Sq.  This  is  the  same  number  as  the  number  of 
available  constants  in  the  asymptotic  solution  at  the  outer  edge  (i.  e. 
and  kn).  Thus  it  will  be  possible  to  match  the  desired  values  C  and  Sq 
at  the' inner  boundary. 


The  value  of  C  must  of  course  be  one  in  all  cases.  The  value 
of  Sq,  which  is  the  entropy  of  the  source  flow,  is  a  matter  of  Choice. 
Since  the  inner  flow  is  isentr.cpic  it  is  useful  to  keep  a  running;  check  on 
the  entropy  (S),  as  given  by 


1 


(4.  20) 


and  the  to^l  energy  H^.  When  these  have  becor^P  nearly  constant  in  the 
process  of  computation,  the  inner  boimdary  has  been,  reached.  Denoting 
the  final  value  of  entropy  by  Sq,  it  can  be  seen  that  it  is  equivalent  to  a 
statement  of  the  chamber  pressure  of  the  Source  flow,  since  from  (4.  20) 
and  the  equation  of  state,  and  Using  the  fact  that  ho  ~  Hq  at  source 
stagnation,  one  has: 


.V% 


>0-1 


^q,  chamber 


There  now  arises  ;Ae  question  of  what  kinds  of  inner  source  flows  the 
solutions  are  capable  of  representing.  It  has  been  noted  that  the  nec¬ 
essity  of  integrating  through  the  dividing  streamline  (stagnation)  may 
place  restrictions  bn  the  choice  of  values  for  k-  and  kg.  in  this 
regard  it  can  be  expected  that  if  the  asymptotic  formulas  ;(^.  19)  are  correct, 
then  any  choice  of  k^,  kg  and  k^  will  yield  a  flow  that  is  regular  at 
stagnation  (i.  e.  thatany  sing^af  behavior  there  is  merely  coznputational 
instability).  This  wOuld  leave  full  freedom  to  match  conditions  on  C  and 
Sq.  The  other  quantities  however.  Such  as  u,  v,  p,  H,  and  S.  cannot 
be  specified,  but  must  be  accepted  as  determined  by  the  process  of  in¬ 
tegration.  It  can  be  shown  for  example  that  if  the  inner  flow  is  to  have 
the  same  entropy  along  streamlines  off  the  axis  (i.  e.  ^  f  0),  then  po 
must  equal  pq;  however  there  are  not  enough  constants  to  enforce  this 
condition.  As  a  practical  inatter  this  lack  is  related  to  the  inconsistency 
in  some  of  the  second  order  relationships  as  discussed  in  the  previous  annual 
report  (ref.  Hill  and  'Habert  (1963)).  The  details  of  the  inner  flow  will  have 
to  await  the  completion  of  a  number  of  solutions  of  the  type  discussed  here. 
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4.  4  Results  of  Numerical  Analysis 


As  noted  in  4.  1,,  preliminary  results  have  been  obtained  using 
an  earlier  method  of  calculation.  Essentially  this  was  to  integrate 
£rom  both  sides  of  the  merged  layer  using  estimated  properties  of  the 
solutions  as  boundary  conditions  and  then  to  attempt  to  match  flow  prop¬ 
erties  at  the  center  (singular  point)  by  adjusting  the  boundary  values. 

The  final  answer  is  obtained  by  applying  successive  corrections  to 
the  eight  undertermined  gradients  and  the  value  of  the  radius  at  the 
inside  boundary.  To  that  end,  it  is  assumed  that  the  value  of  any 
characteristic  quantity,  obtained  when  the  integration  is  stopped  at 
the  singularity  point,  is  a  linear  function  of  the  corresponding  initial 
parameters  with  which  the  integration  has  been  carried  out.  These 
parameters  consist  of  the  four  physically  undertermined  gradients 
at  the  outer  edge,  the  corresponding  four  gradients  and  the  initial  value 
of  the  independent  variable  r  at  the  inner  edge.  At  the  singular  point, 
the  values  of  any  characteristic  quantity  obtained  by  integration  from 
outside  and  inside  is  different.  This  difference  is  a  linear  function  of 
the  nine,  parameters  mentioned  above.  The  new  values  of  the  parameters 
are  obtained  by  computing  the  increment  which  must  be  added  to  each 
in  order  that  the  difference  be  zero.  Since  there  are  nine  parameters, 
a  zero  difference  can  be  imposed  on  nine  characteristic  quantities.  The 
latter  computation  is  carried,  out  by  inverting  a  9  x  9  matrix.  The 
computation  is  repeated  until  the  differences  fall  below  a  srhall  given 
number.  This  scheme  is  sometimes  referred  to  as  the  method  of 
influence  coefficients. 

It  is  assumed  that  at  an  altitude  of  approximately  180  kilometers 
a  missile  is  flying  at  4  kilometers  per  second.  The  Mach. number  is 
6.. 5  and  the  Reynolds  number  is  65.  The  pressure  in  the  combustion 
chamber  of  the  rocket  is  46.  4  atmospheres;  its  temperature  is  3600°K. 
The  exhaust  gases  have  a  ratio  of  specific  heat 'of  1,.  25,  their  molecular 
weight  is  22.  The  result  of  the  computation,  is  presented  on. Figs.  4.2 
and  4.  3.  The  influence  coefficients  scheme  has  been  arranged  in  order 
that  the  followhrg  characteristic  quantities  should  match  at  the  singular 
point: 

p,  V,  u,  c,  h,  ph',  pc',  pv',  pu' 

The  total  enthalpy  matches  by  definition;  the  pressure  term  pg^matches 
too  because 

Pg  =  r  (c)  ph 

The  second  order  term  P2  has  not  been  required  to  match.  It  is  re¬ 
markable  that  the  density,  velocity  components,  enthalpy,  and  concen¬ 
tration  match  with  an  approximation  of  3/  1000  or  smaller  and  the 
gradient  with  an  approximation  of  3  per  cent  or  smaller. 
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VELOCITY  MIXING  RATIO 
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Figure  4.  2.  Stagnation  Point  Merged  Layer  Profiles 
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Figure  4.  3.  Stagnation  Point  Merged  Layer  Profiles 
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DENSITY 


The  sandwich-type  structure  of  the  merged  layer  is  clearly 
illustrated  by  the  results  shown.  There  is  a  thick  outer  shock, 

0.  9<  r<  1.  0  approximately,  in  which  the  thermodynamic  variables  p, 
p  and  h  rise  to  values  very  nearly  equal  to  those  for  a  thin  shock  in 
air  at  M  s  6.  5.  There  is  an  inner  shock  about  half  as  thick,  0.  75<r 
<0.  80,  in  which  the  exhaust  gases  are  compressed.  Here  also  the 
degree  of  corhpression  is  very  nearly  that  which  would  be  computed 
from  the  Rankine-Hugoniot  equations.  The  inner  shock  is  thinner 
because  the  temperature  and  hence  viscosity  of  the  gas  approaching 
it -is  very  much  lower  than  that  of  the  free  stream.  ‘Between  these  two 
shocks  is  a  layer  across  which  the  pressure  is  nearly  constant.  This 
corresponds  to  the  mixing  layer  in  the  high-Reynolds  number  flows 
The  "dividing  streamline,"  v  ■  0,  is  contained  in  this  layer  and 
the  gradient  of  v  is  very  much  smaller  than  in  the  shocks.  At  this 
Reynolds  number,'  however,  the  mixing  extends  beyond  this  region 
into  the  shock  waves  and;  the  variation  of  mixing  ratio  in  the  constant- 
pressure  region  is  only  about  0.  2<c<0.  9- 

The  structure  of  the  merged  layer  at  the  nose  of  a  high-altitude 
rocket  plume  is  now  known  for  one  particular  flight  condition.  It  is 
well  enough  understood  that  qualitative  extrapolation  to  both  higher 
and  lower  altitudes  is  possible. 

We  recomihend  that  additional  and  numerical  integration  of  the 
equations  Be  carried  out  to  cover  the  range  5  <  Re  <  1000.  We  would 
then  be  iii  position  to  continue  the  solution  away  from  the  stagnation 
point  using  approximate  methods.  The  inner  boundary  conditions 
would  be  given  by  the  analysis  of  Section  3. 

The  ultimate  objective  of  this  analysis  is,  of  course,  the  pre¬ 
diction  of  the  observables  generated  by  the  chemical  reactions  be¬ 
tween  the  exhaust  gases  and  the  ambient.  Work  on  this  can  be  started 
immediately  for  those  reactions  for  which  the  chemical  kinetics  have 
been  established. 
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4.  5  List  of  Symbols  Used  in  Section  4 


c 

C 

P  • 
^air 

C 

P 

^gas 

h 


H 

Po’Pz 


r 

R 

Re 


R  . 
air 


gas 


R 


per 


S 

S 

1 

ii 

V 


gas 


V 

Y 

Y 

’gas 

F 

0 

P 

V 
M 


mixing  ratio 

specific  heat  at  constant  pressure  for  air 

specific  heat  at  constant  pressure  for  exhaust  gases 

diffusion  coefficient 
specific  enthalpy 
total  enthalpy 
pressure  terms 

radial  coordinate 

radius  of  the  outside  sphere  (see  Fig.  4.  1.) 

Reynolds  number  based  on  the  radius  of  the  outside  sphere 
and  upstream  conditions 

perfect  gas  constant  for  air 
perfect  gas  constant  for  exhaust  gases 
perfect  gas  constant  for  any  gas 
entropy 

entropy  of  exhaust  gases 

tangential  component  of  velocity 

radial  component  of  ’-elocity- 

specific  volume 

ratio  of  specific  heat 

ratio  of  specific  heat  for  exhaust  gases 

gas  constant  for  the  mixture  of  gases 
angular  coordinate 
viscosity  coefficient 


61 


MC  63-80-RI 


5.  SHOCK  WAVE  FORMATION  BY 
HIGH -ALTITUDE  ROCKET  PLUMES 


5.  1  General 


The  shockwaves  generated  by  a  rocket-powered  missile  in 
high -altitude  flight  are  associated  with  the  plume  rather  than  with 
the  body  of  the  missile.  As  far  as  generating  a  disturbance  in  the 
ambient  is  concerned,  the  rocket  plume  acts  just  like  a  blunt-nosed 
body.  Its  size  becomes  very  large  at  high  altitudes:  typical  diameters 
at  an  altitude  of  200  Km  are  several  kilometers.  Even  though  the 
ambient  mean  free  path  at  these  altitudes  is  very  large  by  laboratory 
standards,,  it  is  still  a  small  fraction  of  the  plume  diameter,  so  that 
continuum  models  of  the  flow  may  be  used,  and  the  formation  of  a 
shock  wave  is  assured. 

As  the  altitude  and  the  size  of  the  mean  free  path  increase, 
the  plume  size  increases  also,  for  a  given  rocket  motor  burning  at 
a  constant  rate.  The  simple  blast-wave  theory  of  Hill  and  Habert 
(1963)  predicts,  however,  that  the  plume  does  not  grow  as  fast  as 
the  mean  free  path,  so  that  eventually  a  free -molecule  model  of  the 
flow  wdll  have  to  be  used  and  no  shock  wave  will  form. 

The  purpose  of  the  work  reported  here  is  the  development  of 
an  approximate  criterion  for  the  upper  limit,  in  altitude,  of  shock 
Svaye  formation.  It  is  generally  agreed  that  this  limit  occurs  in  the 
transition  reginie  between  continuum  and  free  molecule  flows.  This 
regime  is  the  most  difficult  one  to  handle  theoretically  and  as  Talbot 
(1962)  has  pointed  out,  there  exists  no  really  satisfactory  theoretical 
prediction  for  the  limit  of  formation  of  a  well  defined  shock  //ave. 


5.  2  High  Altitude  Flow  Regimes 

Figure  5.  1  copied  from  Levinsky  and  Yoshihara,  (1961)  illus¬ 
trates  schematically  the  flow  regimes  encountered  by  a  re-entering 
hypervelocity  2-foot  diameter  sphere.  The  same  range  of  flow 
regimes  is  encountered,  in  reverse  order,  by  an  ascending  rocket- 
powered  missile.  The  altitude  for  each  regime  is,  of  course,  much 
higher  in  the  latter  case,  because  the  plume  is  much  larger. 

For  the  three  flow  regimes  in  the  lower  row  of  the  figure, 
continuum  theory  applies  and  a  shock  wave  is  formed.  Note  that  the 
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shock  thickness  increases  with  inci’easing  altitud*-  and  dec.^asing 
Reynolds  number.  The  boundary  layer  thickness  <  Isa  incre.ases  and 
in  the  highest.-altitude  continuum  flow  regime  the  '’hock  wave  and 
boundary  layvir  merge  to  gorm  a  single  viscous  layer  between  the  body 
surface  and  the  free  stream.  Hill  and  Habert  (19631  Jhave  discussed 
—the  corresponding  regime  in  the  plume  flow  where  two  shocks  and 
a  mixing  layer  merge  at  the  boundary  between  the  ambient  and  the 
exhaust  gases. 

Consider  now  the  flow  in  the  limit  of  extremely  high  altitudes. 
Here  collisions  of  air  molecules  with  other  air  molecules  are  much 
rarer  than  collisions  between  the  air  molecules  and  the  body  and  may 
be  neglected  in  calculation  of  the  flow  pattern.  This  constitutes  the 
"free  molecular"  regirne  of  flight  in  which  there  i&  no  pile  up  of 
molecules  ahead  of  the  body  and  no  shock  wave  is  formed;  The  corres¬ 
ponding  regime  in  the  plume  flow  is  the  "collisionless  expansion"  of 
exhaust  gases. 

In  between  the  continuum  and  free -molecular  flow  regimes 
lies  a  transition  regime  in  which  collisions  between  air  molecules 
cannot  be  altogether  neglected.  The  molecules  reflected  from  the 
body  cannot  penetrate  indefinitely  far  upstream  without  colliding 
with  other  molecules  in  the  ambient.  As  a  result  of  such  collisions 
there  will  be  a  tendency  for  tpolecules  to  pile  up  in  front  of  the  body 
and  to  form  a  region  of  increased  density.  The  higher  the  ambient 
density,  the  smaller  the  penetration  of  the  reflected  moleculei-  and 
thus  the  thinner  the  high-density  region  in  front  of  the  body.  This 
is  the  process  by  which  a  shock -wave  is  initially  formed  as  a  body  de¬ 
scends  towards  higher  atmospheric  densities.  In  the  case  of  the 
plume  flow  the  transitional  regime  corresponds  to  the  onset  of  colli¬ 
sions  between  the  ambient  molecules  and  those  from  the  rocket  exhaust, 
forming  a  double  shock  containing  a  high-density  region. 

‘  Figure  5.  2  is  a  chart,  again  adressed  to  the  re-entry  problem, 
locating  the  various  flow  regimes  around  a  2 -foot  diameter  sphere 
in  altitude  and  velocity.  Here  Probstein  (1961)  has  broken  some  of 
the  regimes  discussed  above  into  sub-regimes  which  correspond  to 
regions  of  validity  of  various  mathematical  approximations.  Without 
much  explanation  he  locates  the  limit  of  shock  formation  along  the 
line  where  the  shock  thickness  (estimated  at  3  mean  free  paths)  equals 
the  body  radius: 

<Rj^  (5.1) 

5.  3  Criteria  for  Shock  Wave  Formation 


Since  there  is  a  pile-up  of  molecules  in  front  of  a  body  even 
in  free-molecule  flow,  it  is  difficult  to  establish  a  criterion  for 
shockwave  formation  in  terms  of  the  density  gradient  there.  It  seems 
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Figure  5.  2.  Flow  Regimes  for  a  2-foot  Diameter  Sphere  (after  Probstein) 
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preferable  to  examine  this  question  in  terms  of  the  lateral  propa¬ 
gation  of  a  density  perturbation.  Methods  for  calculating  shock  wave 
shape  and  decay  at  very  low  Reynolds  numbers  have  been  discussed 
by  Lighthill  (1956).,  Although  a  specific  calculation  procedure  has 
not  been  set  up  for  either  the  blunt -body  or  the  plume  problem,  Light- 
hill's  results  can  be  used  more  or  less  directly  to  establish  a  criterion. 

The  relevant  conclusion  expressed  in  Lighthill's  article  is  that 
"acoustic  disturbances  decay  according  to  the  power  laws  given  by 
shock  wave  theory  while  their  Reynolds  number  exceeds  about  1,  and 
according  to  the  exponential  laws  given  by  attenuation  theory  when 
IT  is  less,  than  about  -1.  (The  former  decay  laws  are  more  rapid  in  the 
former  range,,  the  latter  in  the  latter. )"  This  means  that  if  R  >  1, 
the  bow  wave  aliead-of  a  body  or  plume  will  propagate  laterally  approxir 
;mately  as  predicted  by  blast-v/ave  theory,  while  if  R  <  1,  the  disturbance 
will  die  out  very  rapi^y  around  the  sides  of  the  body. 

Now  Lighthill's  model  of  thick  shock  waves  with  viscous  dissipation 
is  based  on  Burger's  equation  and  his  Reynolds  numbers,  R,  are  not 
directly  comparable  with  ordinary  Reynolds  numbers  based  on  shock 
wave  thickness.  It  is  convenient,  therefore,  to  consider  the  complete 
N-wave,  consisting  of  the  bow  shock,  the  rear  shock,  and  the  expan¬ 
sion  between  them..  In  this  case  the  Reynolds  number  is  simply  rer 
lated  to  the  ratio  of  shock  thickness  to  the  length  of  the  whole  wave. 

Figure  5.  3  illustrates  schematically  the  N.-wave  generated  by 
a  sphere  in  both  a  high -density  and  a  low-density  flow.  In  the  latter 
the  shock  wave  thickness  is  of  the  order  of  the  sphere  radius.  For 
the  shock  wave  spacing  shown,  this  thickness  corresponds  to  about 
one  quarter  of  the  total  length  of  the  N-wave  and  it  may  be  deduced  from 
Lighthill'  s  calculated  profile  shapes  tha^this  corresponds  to  R  ~2. 

Note,  however,  that  this  argument  is  only  qualitative  since  the  spacing 
between  the  bow  and  rear  shocks  has  only  been  sketched  approximately. 
What  it  does  show  is  that  Probstein's  (1961)  arbitrary  criterion  can  be 
reformulated  in  terms- of  the  initial  Reynolds  number  of  the  N-wave 
generated  by  the  plume  or  :body  and  hence  in  terms  of  the  lateral  rate 
of  decay  of  this  wave.  Note  that  the  value  of  R  in  an  axi symmetric 
flow  decreases  with  increasing  distance  from  the  axis  so  that  a  wave 
which  decays  according  to  the  blast-wave  law  initially  may  eventually 
reach  a  point  where  it  begins  to  decay  much  faster. 

Now  although  the  shock  wave  forms  in  the  transitional  regime, 
estimates  of  its  thickness  can  only  be  obtained  as  extrapolations  from 
the  merged  layer  regime,  where  calculatioyis  are  possible.  These 
calculations  yield  a  thickness  at  the  stagnation  point  and  not  out  where 
the  N-wave  begins  to  form.  Since  the  thickening  of  the  shockwave  be¬ 
tween^  these  two  points  cannot  be  calculated,  we  have  chosen  the  some¬ 
what  conservative  criterion 

Rj^  (5.  2) 
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(a)  HIGH  DENSITY  FLOW 


(b)LOW  DENSITY  FLOW 


Figure  5.  3.  Schematic  N-Waves 
Hypersonic  Sphere 


which,  starting  from  a  completely  different  line  of  reasoning,  only 
differs  from  Probstein  by  a  factor  of  2. 


5.  4  Estimates  of  Shock  Wave  Thickness 

Extensive  calculations  (Liepmann  et  al  (1962))  by  means  of 
the  Boltzmann  equation  of  the  structure  oFplane  shocks  have  show'n 
that  the  Navier-Stokes  model  yields  a  good  approximation  to  the 
thickness,  if  not  to  the  details  of  the  structure.  We  shall,  therefore, 
estimate  shock  thickness  in  blunt-body  and  plurne flows  by  extrapola¬ 
ting  available  Navier -Stokes  calculations. 

Figure  5.  4  represents  the  results  of  numerical  calculations 
of  shock  thickness  in  blunt-body  flows.  Three  points  have  been  taken 
from  Kao's  (1963)  calculations  Of  the  flow  past  a  sphere  and  one  point 
represents  our  own  initial  calculation  of  the  merged  layer  at  the 
nose  of  a  plume.  The  correlation  parameter  is 


M 


oo 


00 

K 


The  inverse  Mach  number  here  accounts  roughly  for  the  variation  be¬ 
tween  restUts  obtained  with  different  values  of  M^.  It  is  consistent 
with  the  idea  that  shock  thickness  should' depend  on  the  mean  free 
path  behind  the  shock  more  strongly  than  on  the  one  in  the  incident 
stream.  ' 


Extrapolation  of  these  data  to  the  critical  value  of  Eq.  5.  2 
predicts  the  criterion 


1 

00* 


^ «  < 
^  - 


0.4 


(5.  3) 


for  the  onset  of  shock  formation. 

The  correlation  parameter  chosen  here  is  essentially  an 
inverse  Reynolds  number  and  it  can  easily  be  shown  that  (5.  3)  is 
equivalent  to 


Re 


Pgo^oo^N 


do 


(5.  4) 


However  the  criterion  in  terms  of  the  mean  free  path  is  easier 
to  use  since  it  is  tabulated  to  higher  altitudes. 
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Figure  5.  4.  Shockwave  Thickness  in  Merged-Layer  Blunt-Bodv 
Flows 
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SHOCKWAVE  FORMATION  CRITERION 


5.  5  Numerical  Results 


The  nose-radius  of  a  high-altitude  rocket  plume  may  be  com 
puted  from  the  blast  wave  formula  of  Hill  and  Habert  (1963) 


R 


(5.5) 


where  D  is  the  drag  of  the  plume  and  Jq  is  a  parameter  in  the  blast 
wave  theory  which  depends  on  Yoo*  Over  the  range  of  Yoo  interest 
the  constant  value  Jq  ■  0.  85  is  a  good  enough  approximation  for  the 
present  purpose. 

The  validity  of  this  theory  in  the  transitional  flow  regime  is, 
of  course,  dubious.  Its  prediction  is  nevertheless  expected  to  be 
approximately  correct.  Neely  (1964)  has  used  it  to  estimate  the 
plume  size  of  the  SYNCOM  apogee  motor  at  an  altitude  of  35,  000  Km 
and  found  agreement  to  within  half  an  order  of  magnitude  with  Baker- 
Nunn  photographs. 

The  shock  thickness  correlation  parameter  may  therefore  be 

written 


Figure  5.  5  illustrates  the  variation  of  this  parameter  with  altitude  for 
two  values  of  plume  drag.  Note,  that  the  plume  drag  value  D  *  12,  800 
lbs  corresponds  to  the  Atlas  or  Titan  sustainer  engine.  Also  shown 
in  Fig.  5.  5  is  the  limiting  value  of  this  parameter  for  shock  wave 
formation,  criterion  (5.  3).  It  predict.*^  shock  wave  formation  to  very 
high  altitude,  above  400  Km  for  the  Atlas  and  Titan  engine. 

The  continuum -flow  criterion  used  by  Hill  and  Habert  (1963) 

was 


This  can  be  located  in  Fig.  5.  5  only  if  a  value  of  Mo©  chosen.  It 
has  been  drawn  for  Mao  *  6,  which  is  typical  in  the  altitude  range  of 
interest.  Note  that  it  lies  below  the  shock  wave  criterion  by  about 


70 


MG  63-80-Rl 


71 


MC  63-80-Rl 


one  order  of  magnitude.  According  to  Probstcin's  (1961)  chart  illus 
trated  in  Fig.  5.  2,  this  spread  is  too  ilarge;  a  less  severe  criterion 
for  continuum  flow  would  be  more  appropriate. 

5.  6  Conclusions 


The  development  in  this  section  of  a  criterion  for  shock  wave 
formation  by  high>altitude  rocket  plumes  is  based  on  estimates  of 
the  Reynolds  number  of  the  N-wave  formed  by  the  bow  and  rear  shocks 
generated  by  the  plume.  The  following  conclusions  may  be  drawn: 

1.  A  more  precise  expression  of  the  criterion  developed 
will  require  a  better  model  of  the  rear  shock  in  a  plume 
flow. 

2.  Given  the  length  of  the  N>wave,  the  new  Criterion  may 
easily  be  related  to  the  more  fainiliar  criteria  given  in 
the  literature  in  terms,^Qf  the  ratio  of  shock  thickness 
to  nose  radius, 

3.  Extrapolation  of  Navier-Stokes  solutions:  for  both  sphere 

and  plume  flows  has  provided  simple  estimates  of 
6g/Riif  in  terms  of  the  parameter  which  is  of 

the  form  of  a  reciprocal  Reynolds  number. 

4.  Using  plume*size  estimates  given  by  blast-wave  theory 
it  has  been  shown  that  a  shpek  wave  will  be  fornied 
throughout  the  burning  period  of  typical  ICBM  engines. 
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